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ABSTRACT 

Radially inhomogeneous gamma-ray burst (GRB) jets release variable photosphcric 
emission and can have internal shocks occurring above the photosphere. We generically 
formulate a photospheric emission model of GRBs including Compton up-scattered 
photospheric (UP) emission off the electrons (and positrons) in the internal shocks, 
and find that the photospheric emission may correspond to the traditional (Band) 
component at < 1 MeV and the UP emission to the high-energy emission observed 
by Fermi/LAT for some GRBs at > fO MeV. The two components can be separate 
in the spectrum in some cases or can mimic a smooth broad Band spectrum in other 
cases. We apply our formulation to the well-studied long and short LAT GRBs, GRB 
080916C, GRB 090902B, and GRB 090510, and typically find reasonable parameters 
for fitting the time-binned spectra, although fine tuning of several parameters is re- 
quired. The observed delays of the high-energy emission with respect to the MeV 
emission which are large compared to the variability times are unlikely to be due to 
simple kinematic effects of a non-evolving jet. These delays may instead be attributed 
to the temporal evolution of the physical parameters of the jet, and thus the delay 
timescales could provide a potential tool for investigating the structures of GRB jets 
themselves and their progenitors. The difference of the delay timescales of long and 
short GRBs inferred from the Fermi data might be due to the differences in the pro- 
genitors of long and short GRBs. Some other properties and consequences of this 
model are discussed, including temporal correlations among the prompt optical, the 
soft X-ray, and the distinct high-energy component as well as the Band component. 
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1 INTRODUCTION 

The origin of gamma-ray burst (GRB) prompt emission re- 
mains unclear. It has been mainly observed in the energy 
range of 10 keV — 1 MeV, and most of the spectra are fit- 
ted by a simple broken power-law functi on (so-called Band 
function) or a cutoff power-law function |Preece et al.ll2000l ; 
iGhirlanda et alj|2002l : iKaneko et al.ll2006r ), the light curves 
being typically very variable. Many models have been pro- 
posed to explain these properties. Most of them invoke a 
relativistic jet that is energized by a newly born compact 
object. Currently three types of emission mechanisms in the 
relativistic jet are being actively discussed; photospheric, 
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leptonic synchrotron, and hadronic emission models. The 
first models assume that the thermal energy stored in the 
jet can be radi ated as prompt emission at the Thomson pho- 
tosphere (e.g., Paczynski 1986; G oodman! Il986l ; iThompsonl 
1994; Mes zaros fc Reesll2000l ). Here the thermal energy in 
the jet can be produced by the particle and/or magnetic en- 
ergy dissipation between the explosion center and the photo- 
sphere. The second and third models assume the Thomson- 
thin region as the prompt emission site. In the second 
models, the electrons (and positrons) accelerated by the 
shock dissipation of the kinetic energy or by the magnetic 
energy dissipation radiate synchro tron and synchrotron- 



self-Compton (SSC) em i ssions ( e.g., Rees fc Meszarod 



Daigne k, Mochkovitchl 1 19981: iKumar fc McMahonl 



1994; 



2008; 



Spruit et al.ll200ll : lLvutikovll2006r i, while in the third mod- 
els the accelerated protons and secondary particles induced 
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by the photopion cascade process produce synchrotron 
and inverse Compton (IC) emissions fe.g.. IVietril 1 19971 ; 



iBottcher fe Dermeri 1 19981 : lAsano et all l2009bf l. Clarifying 
the prompt emission mechanism and the physical quanti- 
ties at the emission site would help understand the nature 
of the relativistic jet and the compact object that energizes 
the jet. For ot her mo d els an d more gen e ral di scussio ns, see 
recent reviews iPiranl (|2004h , iMeszaros! (|2006l l. and IZhangj 
l|2007l ). 

GRBs were only sparsely observed in the > 10 MeV 
energy range, until the Fermi satellite was launched 2008 
June 11. Now the GBM (8 keV - 40 MeV) and the LAT 
(~ 20 MeV — 300 GeV) detectors onboard Fermi provide 
extremely broad energy coverage with good temporal res- 
olution for GRBs. The Fermi observations will put further 
constraints on the above three types of prompt emission 
models. During its first 1.5 yr routine operation, t he LAT 
has d etected 14 GRBs. Those are summarized in iGranotl 
(|2010h . Most of their spectra are fitted by a Band function 
even up to ~ 10 GeV, while at least 3 GRBs (GRB 090510, 
GRB 090902B, and GRB 090926A) have additional distinct 
spectral component at > 10 MeV. Those additional compo- 
nents are fitted by single power-law function. These 3 GRBs 
are among the brightest GRBs detected by Fermi. This sug- 
gests that such a high-energy component may be very com- 
mon and we can clearly detect such a distinct high -energy 
component only in bright LAT GRBs l|Granotll2010r i. 

Fermi also revealed that the high-energy emission (> 
100 MeV) of most LAT GRBs is delayed behind the on- 
set of the MeV emission, and the high-energy emission of 
many LAT GRBs lasts longer than the MeV emission, show- 
ing power-law decays, which are typically detected until 
~ 10 2 — 10 3 s. The delay times in the cosmological rest 
frame are ~ 1 s for long GRBs and < 0.1 s for short 
bursts GRB 081024B and GRB 090510. Although some de- 
lays of the high-energy photons are just caused by the flux 
increases above the LAT detection threshold without a spec- 
tral change, others must clearly be attributed to the spec- 
tral changes of the Band component and/or the onset of 
the distinct spectral component. We note that the distinct 
components of the 3 GRBs are delayed. 

The origins of the distinct spectral components and the 
onset delays in the high-energy range have been actively de- 
bated. It has been proposed that the high-energy emission 
can be attributed to the external shock, which is made by 
the interaction of the jet with the ambient medium, and 
the onset d elays can correspond to the times for the jet de- 
celerations (Kumar & Barniol Duran 2009; Ghisellini ct al. 
2010T) (see also iGranot fc G~u ctta 2003]; iPe'er fc Waxmanl 



2004). However, the observed high-energy emission usually 



has a strong variability and it often correlates with the 
MeV emission, which is at odds with an external shock 
origin. Especially for GRB 090510, whose long-lived high- 
energy emission can be explained by the external shock 
synchrotron emission together with the optical and X-ray 
afte rglows detected from ~ 100 s after the burst trig- 
ger (|Pe Pasquale et al.ll2010l ; lKumar fc Barniol Duranl2010l; 
ICorsi et alj|2010l ). it is found that the high-energy emission 
in the prompt phase is muc h brighter than that produced 
by the same external shock (|He et al.ll2010l ). A similar con- 
clusio n has been obtained for GRB 090902B i|Liu fc Wand 
2010). As for the internal shock models, it is not simple 



to explain the LAT onset delays which are larger than 
the variability times cales (t v ) by the de layed brightening 
of the SSC emission (lAbd o et a l. 2009b , cj; | Ackermann et al.l 
l2010al ; ICorsi et al.ll2010l ; iDaigne et al.M201ll ). The spectral 
index of the distinct high-energy component different from 
the Band low-energy spectral index, seen in e.g., GRB 
090902B, may not be explained simply either in these mod- 
els. To overcome these problems, we have proposed an ex- 
ternal inverse Compton (EIC) component (i.e., the emis- 
sion produced by up-scattering the photons incident from 
outside the shock) in addition to the sync hrotron and 
SSC c omponents in the internal shock model l|Toma et al.l 
l2009bh ; this model, however, requires an extreme value 
of the microphysical parameter e_g ~ 10 -5 , which does 
not seem common. Other alternative mechanisms such as 
hadr onic emission mechanisms require huge isotropic ener- 
gies (lAsano et a l. 2009a; Wa ng et alj 120091 ; iRazzaque et alj 
l2010h . while magnetically-dominated jets have not been 
explored explicitly enough for detail ed comparisons with 
Fermi/LAT observational res ults (see IZhang fc Pe'erl 120091 ; 
lFanll2010l ; IZhang fc Yanll201 lT ). 

In this paper we concentrate on a photospheric emis- 
sion model. Such models have been briefly discussed for in- 
terpreting the relation of the MeV emission with the high- 
energy emission. In this scenario, the photosphere produces 
a Band emission component which peaks around ~ 1 MeV, 
and this cannot generally produce the high-energy emis- 
sion because of large opacity for e ± pair creation. The 
high-energy emission in these models ma y instead arise in 
a dissipation region at a larger r adius (|Gao et al.1 120091 ; 
lBeloborodovll2010l ; iRvde et alBoioh . In this paper, we sub- 
stantially extend this idea and discuss its consequences in 
greater depth. We focus on the case in which the energy of 
the jet is mainly carried by photons and baryons (where 
magnetic field energy is subdominant), and the dissipa- 
tion mechanism at large radius is internal shock. We show 
that in typical cases the photospheric emission is efficiently 
Compton scattered by the electrons in internal shocks out- 
side the photosphere, and find that the up-scattered photo- 
spheric (UP) emission is a good candidate for the observed 
high-energy emission in the prompt phase. This Comp- 
ton scattering is a type of the EIC scattering, which is 
thought to be important for the high-energy emission of 
blazars (e.g..lSikola et a} . 1994; Do rmer fc Sc hlickciser1 [l993l ; 
iBrune tti 2000; Ghisellini fc Tavecc hio 2009), and could also 
be important in the internal and external shocks in GRBs 
teeloborodoyl2005l;IWang fc Mesza^l200d ; lFan et alj2008l : 
iToma et al.1 l2009bl ; iMurase et all l2010al lbV We make a 
generic formulation of the spectral types of the radiation 
from the photosphere and the internal shock including syn- 
chrotron and SSC emission for the cases of the efficient EIC 
scattering, and clarify necessary conditions for the photo- 
spheric and UP components to be dominant in the energy 
range of 10 keV — 10 GeV, rather than synchrotron and 
SSC components. This formulation is applied for the well- 
observed LAT long and short GRBs, GRB 080916C, GRB 
090902B, and GRB 090510, and we typically find reason- 
able parameter sets for which the data can be explained by 
the photospheric and UP emission. Our model fits indicate 
that the observed delayed onset of the LAT emission may 
be interpreted as the parameter evolution of the GRB jet. 

The UP emission ceases at the end of the prompt MeV 
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(photospheric) emission, and the subsequent long-lived high- 
energy emijjsionjma^^ 

(e.g., iDe Pasauale et all l20ld : iKumar fc Barniol Duranl 
I201CJ : iGhisellini et alj l20ld )~~A different origin for the 
prompt and late LAT emission is in fact implied by the 
change in its beha vior across this transition l|He et al.ll2010l : 
iLiu fe Wan3l201Ch . In what follows, we concentrate on the 
prompt emission phase. We discuss the general temporal 
and spectral properties of the radiation from the photo- 
sphere and the internal shock of the GRB jet in Sections 
HO and H and perform the case studies of GRB 080916C, 
GRB 090902B, and GRB 090510 in Section [5] Some of the 
implications of our results are discussed in Section [6] 



2 PHOTOSPHERIC EMISSION 

We consider that the jet is accelerated by the thermal pres- 
sure, and the magnetic field energy is subdominant. We de- 
rive the luminosity and temperature of the photospheric 
emission from the jet and the remaining kinetic luminos- 
ity of the jet above the photosphere, which is the luminos- 
ity budget for internal shocks , according to the standard 
fireball model (|Paczvnskilll986l. [l990l; IShemi fc Piranlll99d; 



Meszaros et al.l 19931 ; iMeszaros fc Reesl l200dT iNakar et al.l 
20051 ). We assume that the jet material is optically thick 



and the energy is dominated by thermal radiation at the 
base, r — r a , from where it expands in the roughly adia- 
batic condition, i.e., without strong conversion of the kinetic 
energy into the thermal energy affecting the adiabatic ex- 
pan sion dynamics (fo r a co ntinuous strong dissipation case, 
see iRees fc Meszarosl |2005h . The jet material may have a 
Lorentz factor r a 1 at the base. In particular, if the GRB 
originates from a massive stellar collapse and the jet suffers 
strong dissipation of its kinetic energy due to interaction 
with the stellar envelope, the front portion of the jet may 
have r a ~ 10 10 - 10 11 c m and Y a ~ 10 - 10 2 (see recent 
nume rical simulations bv lMorsonv et al1l2007l : lLazzati et al.l 
2009, and discussion in Section^. In the absence of a strong 
dissipation the base radius may be as small as a few of 
Schwarzschild radius of the central compact object, so that 
r a ~ W ( '(A4 c /Mq) cm and T a ~ 1, where M c is the mass of 
the central object. 

The observer-frame temperature at the base is 



T a = 



( LYl 

\Airrlci 



1/4 



-1/2 



MeV/fc, 



(1) 



where L is the isotropic equivalent luminosity of the jet, 
and a ~ 7.56 x 10~ 15 erg cm" 3 K~ 4 is the Stefan con- 
stantQ For general cases of GRBs, the adiabatic index of 
the jet material below the photosphere is found to be 4/3 
(see below). The dynamics of the expanding jet can be sum- 
marized as r oc r and T = T a for r ^ r s , and V = r\ and 
T oc r~ 2 / 3 for r ^ r s , where the Lorentz factor satura- 
tion radius r s = r a rj/Y a - Here we have defined r) = L/Mc 2 , 
where M is the isotropic equivalent mass ejection rate. We 
see that the ratio of the radiation and baryon entropy den- 
sities ~ aT' 4 /n'kT' is constant throughout the evolution 



1 We adopt the notation Q x = Q/10 x in cgs units throughout 
this paper. 



below the photosphere and is much larger than unity if 
r, > kT a /{m p c 2 ) ~ 2 x 10" 3 ^Va,7/r a )- 1/2 (where ri 
and T' is the baryon density and fluid temperature in the 
comoving frame). Thus the assumption of the adiabatic in- 
dex = 4/3 is validated. 

The e ± pairs drop out of equilibrium at T' ~ 20 keV, at 
a radius r~r p = lx 10 9 Ly^{r a ^ /Ta) 1 ^ 2 cm. We consider 
the case in which the jet carries enough baryons to provide 
an electron scattering photosphere above r p (i.e., the case of 
rj < 5 x 10 B Ll^ira^/Ya) 112 ). Then the photosphere radius 
r p h is defined as r — OT"n'r/{2T) = 1, where ot is the 
Thomson cross section and the comoving electron density is 
given by n' = L/(47rr 2 m p c 3 7?r), if there is no ad ditional e 
pair creation in the jet (lAbramowicz et al1ll99ll ). 

If no energy dissipation takes place around the photo- 
sphere, the photons and particles are fully thermalized, typi- 
cally with the comoving temperature T ph ~ 0.1 — 1 keV, and 
then the emerging radiation from the photosphere is a black- 
body with the observer-frame temperature T p h — r p hT ph 
(where T p h is the Lorentz factor at r = r p h). However, some 
energy dissipation processes are expected to occur around 
the photosphere, which can make the emerging radiation 
spectrum deviate from a blackbody. Such processes have 
been demonstrated convincingly by recent analytical and 
numerical work: Both below and above the photosphere, a 
fraction of the electrons can acquire high-temperature distri- 
bution (3> Tp h ), either due to inter nal shocks and/or interac- 
tion of a jet with stellar envelope dEichler fc Levinsonll2000l; 
Rees fc Meszarosl [20051; iPe'er et all 120051 ; iThompson et all 
20071 ; lLazzati fc Bcgclm arJ I2010T ). dissipation of excited 
plasma waves (Io ka et al ] l2007l ^_ or nuclear collisions be- 
tween protons and neutrons 1 Beloborodovl I2010T I. or in 
magnetically-domin ated jets due to scattering with turbu- 
lent Alfven waves dThompsonlll994|) or heating caused by 
magnetic reconnection (iGianniosfe ooe). Multiple IC scatter- 
ings of the thermal photons by mildly relativistic electrons 
(with 7 ~ 1) can create a power-law tail extending from 
the thermal peak, which saturates at the comoving-frame 
energy e' ~ m e c 2 as the Klein-Nishina limit or by the direct 
Compton cooling, which is boosted to e ~ 0.1 — 1 GeV in 
the observer frame. The dissipation processes also can pro- 
duce a relativistic electron population, which IC-scatter the 
photons to even higher energies, the limitation being the 77 
absorption against lower energy photons. The low-energy 
spectral slope of the emerging photons should not be sig- 
nificantly ^_clum|2jd_2nsuch^ 

a ph =± 1 (|Beloborodovll201ol ; lLazzati fc Begelmanll201Ch . 



In our baryonic jet model, we assume that the dis- 
sipation of the kinetic energy is not so strong that the 
jet dynamics is described by the adiabatic evolution, and 
the temperature of the emerging photons is given by T p h. 
This implies that the energy flux in the possible non- 
thermal part above e p h should be smaller than that in 
the main thermal part, leading to an upper limit on the 
photon index /3 p h of the possible power-law tail. Let us 
simply write down the spectral shapes of the pure ther- 
mal component and the thermal plus a non-thermal tail 
spectrum as F e th = 0.84(e/fcT ph ) 3 /[exp(e/fcT ph ) - 1] and 
F e nth = (e/e ph ) 2 min[l, (e/eph)' 3 ? 11 " 1 ], respectively, where 



I ph 



£ p h/4, and these are normalized as F E 



F" 
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at s = s ph . El Then we have F £ th de = 0.84(7r 4 /60)e ph 
and F? th de = [(1/3) - (/3 ph + 2)- 1 ]e ph . The condition 
Jo°° F e nth de < 2 x J 00 F e th de provides a constraint /3 ph < 
—2.5. The value of /3 p h is expected to be determined not 
only by the ratio of the photon to kinetic energies but also 
on the electron energy distribution which depends on the 
type of a dominant dissipation process, which we may not 
specify at the current stage. Hereafter we treat /3 p h as a free 
parameter ranging within the above constraint. 

In such dissipation processes copious e + e~ pairs may 
be created, which increase the Thomson photosphere ra- 
dius. We parameterize the number density of electrons plus 
positrons (which we will call 'leptons' hereafter) near the 
photosphere as n[ = IZn' , where 1Z ^ 1. Typically we 
can have 1Z ~ 10 — 10 2 , in which case the inertia is still 
dominated by baryons. Thus the photosphere is defined as 
n = o T n\r/{2T) = 1. 

We have r p h < r s in the low baryon load case r\ > 77*, 
where 



Tf» 



8nr a m p c 3 



1/4 



2-8 xlO 3 L^(^)~ 17 V 4 .(2) 



In this case, most of the luminosity is radiated at the 
photosphere. Leptons can frequently interact with pho- 
tons even above the photosphere, and this Compton drag 
acceleration determin es the remaining kinetic luminosity 
ijMeszaros et al.lll993l ). The emission radius, luminosity, and 
peak energy of the photospheric emission and the final bulk 
Lorentz factor and remaining kinet ic luminosity of the jet 
above the photosphere are given by ijMeszaros fc ; Reejl2000l ) 



-1/3 



V* / V 



Lph. ~ L, £ p h — 4 kT a , 



(3) 



On the other hand, in the high baryon load case r\ < 77*, 
we have r p h > r s . In this case most of the initial thermal 
energy has been converted to the kinetic energy below the 
photosphere. The corresponding quantities are given bj(f| 



77* / 77 



Lph — L \ — 



8/3 



E ph ~ 4 kT a 



■n* 



8/3 



, r f = 77, L k » L. 



(4) 



For typical GRB jets with parameters L ~ 10 — 
10 54 erg s _1 , the photospheric emission can be bright with 
e p h ~ 100 keV — 10 MeV since the jets can have 10 _1 < 
r a ,7/Y a % 10 2 and 77 > 77, in principle. 



2 More accurately, one should use = 
0.7(e/fcT ph ) 3 /[exp(£//cT ph ) - 1] where T ph = e ph /3 for the 
pure thermal component since e = 4fcT p j 1 is the peak energy 
of the sF e spectrum instead of the F s spectrum. However, we 
obtain the same constraint /3 p h < —2.5 also from this calculation. 

3 We do not discuss the case r p h > r ap = 2Wr/ 2 , for which the 
radial spreading of shells (with the initial wid th W) is significant , 
and the scaling law is different, r p ^ ex. 77 -1 / 2 jNakar et alj|200sh . 



3 TEMPORAL PROPERTIES 

The GRB jet may be thought to consist of many succes- 
sive shells with initial radial widths W <§C r p h. This implies 
that the photospheric emission is temporally variable. After 
emerging from the photospheric regions, collisions between 
the ejected shells can produce internal shocks at n 3> r p h. 
As we will show below, the leptons in the internal shock 
of two given shells can up-scatter their own photospheric 
emission, which was produced as they emerged from the 
photosphere. The UP emission can appear as a high-energy 
spectral component in the LAT energy range (see Section [4}. 
The temporal properties of the emission in this model can 
be exemplified using a simple two-shell collision model (see 
Figure [1]). For simplicity, we only consider the cases in which 
the two shells are both in the regime 77 < 77* or both in the 
regime 77 > 77,. 

The duration of the observed photospheric emission 
from one shell is given by the light crossing time of the shell 
plus the angular spreading time, 



W 

Stph ~ h 



7-ph 

2cr 2 ' 



(5) 



This can be rewritten as 8t p h ~ (W/c) + r p h/(2cr] 2 ) ~ 
(W/c) + (r ph /ri)t v for 77 < 77,, and St ph ~ (W/c) + 
r ph /(2cF 2 h ) ~ (W/c) + (77/V t )(77/77,) 1/3 ^ for 77 > 77,, where 
r p h = 77»(77/t7*) _1/ ' 3 , and 77 = r a r). t /Ya is the radius at 
which the Compton drag acceleration ceases. 

Now we consider a two-shell collision. We assume that 
the two shells have similar initial radial width W , and the 
initial separation is ct v . First consider the high baryon load 
case 77 < 77*. The rapid shell with Lorentz factor T r — r\ r 
catches up with the slower shell with Lorentz factor T s = ?7 S 
(Hereafter subscripts V and 's' denote the quantities of the 
rapid and slow shells, respectively). Setting the zero time 
in the lab frame as the time when the rapid shell is at r — 
r p h, the collision time is given by t; = ct v /[c(/3 r — /3 S )] ~ 
2ct v rj 2 , where we have used an approximation 77 3> 1. The 
dissipation of the total kinetic energy is most efficient whe n 
the masses of the two shells are similar (c.f.. jPiranl 12004 ) . 
For typical cases, the ratio of the two Lorentz factors may 
be a_L = ri r /rj a < 5. Then the dissipation efficiency is given 
by e d < 0.25. 

The photons which arrive at r = n at time t^ti can be 
up-scattered by the energetic leptons in the internal shock. 
Let us define a length l ns , measured from the front end of the 
rapid shell, as the radial scale of the region from which the 
photospheric photons do not interact with the leptons in the 
internal shock. Then we have r p h— ins+ctj =Ti = r p h+cP r ti, 
which gives 



l ns = c(l - p r )ti 



ct v . (77 < 77, ) 



(6) 



We find that inverse Compton scattering is efficient for the 
case l ns < W/2 and inefficient for the case l ns > W/2. 

A similar argument can be made also for the low 
baryon load case 77 > 77,. In this case the shells are ac- 
celerated similarly and cannot collide with each other be- 
low rf = r p h (77/77, ) 1//3 . The internal shock radius is given 
by Ti w 2ct v Ti 2 s , and we consider a case n 3> r/, r . We have 

lns = Cti-Ja* Cfirdt W (Tf^ a /rf^ r )ct v + (r p h, r T p ^ r ~ rf, r T)~r) ■ 

The second term represents a correction due to the velocity 
of the shell lower than 77, , r at r < r/, r . The second term is 
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Figure 1. Interaction of the photospheric emission from two relativistically moving shells with the internal shock that the given shells 
themselves give rise to. The photons from the shaded region can be scattered by the leptons in the internal shock. 



rewritten as r/,r?7*~r[( ? ?r/?7* 1 r) — l]i which cannot be dom- 
inant unless 77 2> r\, , so that l ns can typically be written as 
Equation ((6| also in this case. 

We have two parameter regimes for the scattering effi- 
ciency: 

(i) Efficient scattering regime, l ns < W/2. This includes 
the typical case W ~ ct v . In this regime we may observe 
bright UP emission. The duration of the photospheric emis- 
sion can be written as 5t p h ~ W/c, since we have assumed 
that n ^> r ph for r] < 77, and n ^> rj for r\ > rj t . The onset 
of the UP emission is correlated with that of the photo- 
spheric emission pulse released from the rapid shell and de- 
layed from that of the photospheric emission pulse released 
from the slow shell by ti ag = (W + ct v + l nB )/c ~ (W/c) + t v . 
The duration of the UP emission is given by the duration of 
the seed photons plus the angular spreading time 



5tu 



W-lns 



+ 



3ri 
2cTl 



W _i 
— + 3a L t v - 



(7) 



The factor of ~ 3 in the second term means that the 
UP emission has an anisotropic energy distribution in the 
comoving frame, being brightest at an angle 9 ~ 1/T m 
from the line of sight in the observer frame ( c .f.. iBrunetti 



2000; 



Wang fc Meszarosll2006l ; iFan et al.|[200a : iToma et al 



200901 ). where P m ~ (r r r s ) 1/2 is the Lorentz factor of the 
merged shell. 

(ii) Inefficient scattering regime, l ns > W/2. In this 
regime, typically we may not observe bright UP emission. 
A caveat is that this regime includes a case l ns ~ W + ct v , 
in which we have a third shell behind the rapid shell of the 
two given shells apart by ct v , and the photospheric emission 
from the third shell can be scattered by the leptons in the 
internal shock of the two given shells. This condition reduces 
to a~[ 2 ct v ~ W + ctv, which is satisfied when t v <C t v . 

If the GRB jet is in the efficient scattering regime right 
from the start, when it first emerges, the delay timescale 
between the onsets of the first photospheric emission (in 
the MeV energy range) and the first UP emission (in the 
high-energy range) is ~ ti ag ~ (W/c) + t v , which is compa- 
rable to <5t p h or tv However, Fermi observations show that 



the LAT onset delays are much larger than the variability 
timescale apparent in the MeV energy light curves both for 
long and short GRBs. Therefore it is unlikely that the large 
delays of the high-energy emission onsets observed in many 
GRBs are due to the above simple kinematic effect of the jet 
whose physical parameters do not evolve. In our model, the 
large delays may be interpreted as the timescale on which 
the physical parameters of the jet temporally change, e.g., 
W and t v change from the inefficient scattering regime into 
the efficient scattering regime, or L, r a /T a , r\, and 1Z change 
from the regime 77 > 77* into the regime 77 < 77, (see the fol- 
lowing sections for details of the latter possibility). In these 
cases the UP emission can start, being still dim compared 
with the photospheric emission, at a time ~ ti ag after the 
first photospheric emission, and can become bright with the 
delay as observed. 



4 SPECTRAL PROPERTIES 

As shown above, the photospheric emission of the shells 
can be efficiently up-scattered by the leptons in the inter- 
nal shocks of the same shells. Here we derive the generic 
broadband spectrum of the observed radiation arising in the 
internal shock including synchrotron and SSC emission for 
the efficient scattering regime. 



4.1 Case of 77 < 77, 

In this case the kinetic luminosity that can be dissipated 
into radiation by the internal shock is ~ L. At r > 
r p h, we have arn' l r/(2rj) < 1, and the pair annihilation 
timescale is longer than the expansion timescale. Thus 
the pair population freezes out and we can write n\ = 
TZL/(4iYr 2 m p c 3 r/ 2 ) even well above the photosphere. Here- 
after we focus on the photospheric emission from the rapid 
shell, which is up-scattered by the leptons in the internal 
shock, so that 77 denotes 77,. (The photospheric emission 
from the slower shell may be much dimmer). The internal 
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shock radius is estimated as r, ~ 2ct v r) 2 = 2ct v rfa~^ ~ 
2 x 10 13 »73t t ,,_2(aL/5) _2 cm. Our assumption rj 3> r p h re- 
duces to 



i)»5x!0 2 Lgf ^ 1 /5 t 



-1/5 

u '- 2 V 5 J 



2/5 



(8) 



We assume that all the leptons participate in the non- 
thermal acceleration process in the internal shock and have 
the injected spectrum dn' l i /d'y — C*y~ p for 7 ^ 7™, where 
C is a constant. Then we have C7m _p /(P — 1) = nj(r = r^) 
and Cm e c 2 7j l ~ ! '/(p — 2) = Ltd^e / (4:irr 2 crj 2 ) , where e e is the 
fraction of the dissipation energy that is carried by leptons 
in the internal shock. (We consider a relatively weak internal 
shock, i.e., at < 5, where we can neglect the change of the 
Lorentz factor of the jet for simplicity.) These lead to 



= 7ft e<*ee - 4 ^! — - ,f(p), 

m e p - 1 V 0.1 y 



(9) 



where /(p) = 13(p— 2)/[3(p— 1)]. If only a fraction of leptons 
are accelerated, 7 m is larger than this value. 

We also assume that a fraction es of the dissi- 
pated energy is carried in the form of magnetic fields, 
U' B = B' 2 /(8tv) — LtdtB 1 ' (^nr 2 cr/ 2 ) . The characteristic syn- 
chrotron energy of leptons with Lorentz factor 7 is given 
by 



yn(7) 



3heB' 2 

7 7) 

4nm e c 



Cv1 „-1 2 r l/2 -2,-1 / a L \ 2 I tdtB 

6 x 10 7 L B ' 3 % t^.a ^— J 



eV(10) 



The leptons emit the UP emission, synchrotron emis- 
sion, and SSC emission. The comoving energy density of the 
synchrotron emission can be written by 



U' 

w syn 



t'dyn-^VTCU' B 



2dn[ 1 
7 7j7- d 7 = xU B , 



(11) 



where t' dyn ~ n/(2cq) is the dynamical timescale of the in- 
ternal shock and dnj/efry is the lepton energy distribution 
averaged over the dynamical timescale. The quantity x is 
calculated as 



x=-a — [ 2 dn ' 1 tfrv - ~ ^ 
3 T 277 J 7 d 7 3(p-2) 



ri,i-f m -ych(j m ,jc), (12) 



where 7 C is the cooling Lorentz factor of the leptons and n t i 
is the Thomson optical depth at n (which can be written as 
T i,i — fphAi in the case of r\ < 77*). We have assumed p > 2 
and defined the function /i(7m,7c) as h — 1 for 7 C <C 7m, 
/i = p/(p — 1) for 7 C w 7™, h = (7c/7m) 2 ~ p /(3 — p) for 
7m < 7c an d P < 3, and h = (p - 2)7 m /[(p - 3)y c ] for 
7m <C 7c and p ^ 3. Similarly we have the SSC emission 
energy density U' ssc = a;(7 s ' yn , where we have assumed that 
the Klein-Nishina (KN) effect is not significant. This is valid 
for the parameters adopted in the case study of observed 
LAT GRBs in Section [5] We also have the UP emission 
energy density U^ p = xU ph , where we have assumed that the 
KN effect is not significant also for the UP emission. This 
is found to be valid for observed LAT GRBs in Section [S] 
We also have neglected the anisotropy of the UP photon 
energy distribution, for simplicity. This anisotropy leads to 
the reduction of the observed UP luminosity averaged over 
a pulse than the isotropic assumpt ion by a factor of ~ 2 
jFan et alj|2008l : iToma et alj|2009rj ). which we will neglect 



below. By using the relation, e.g., L p h = 4Tvr 2 cr/ 2 U ph , we 
have 

L Byn ~ € d e B xL ~ kxL pb , (13) 
L BBC ~ dd^BX^L ~ kx 2 L p h, 
where we have defined 

k 'm^' « 141 

The photospheric luminosity is given by L p h ~ L(rj/r/ t ) a ^ 3 
(see Eq. 0} . 

We require to estimate the cooling Lorentz factor 7 C for 
specifying the lepton energy distribution averaged over the 
dynamical timescale and the various emission luminosities. 
Since the cooling rate for one lepton with Lorentz factor 7 
is P(j) = (4/3)a T cy 2 (U' B + U Byn + U ph ), the cooling Lorentz 
factor can be estimated by 7cm e c 2 = P(y c )n/(2cr)). This 
reduces to 

3m e 7^ 1 , s 

7c ~ 4m p r i , i ( ?? / J ? < ,)8/3 jfe(i + x ) + 1 ' [ > 

If k > 1, we can take k(l + x) + l « k(l + x). Then Eq. (fT2)l 
reduces to x = (e e h/ es) / (1 + x), and we have 



e e h/eB, 



(x> 1), 
(x< 1). 



(k » 1) 



(16) 



This is a general result for the case in which the radiation is 
do minated by s y nchro tron and SSC emissions, as developed 
by ISari fc Esinl (|200ll 'l. In the first case, x 3> 1, we have 
kx 2 ^> kx ^> x ^> 1, so that the order of the four emission 
luminosities is found to be L ssc 3> L syn 3> L up ^> L p h- The 
condition k 2> 1 and x 3> 1 is equivalent to ed^ e h 2> edts 2> 
(v/v*) S ■ The second case, 1 < 1, should be divided into 
three sub-cases; (kx 2> 1 and kx 2 3> 1), (kx 3> 1 and kx 2 <C 
1), and (kx <C 1). The orders of the luminosities in these 
cases are summarized in Table [1] which are labeled as cases 
5, 6, 7, and 8, respectively. In the Table we define Q = 
(?j/?7*) 8 ^ 3 , £ = edteh, and B = td^B for clarification. 

If k < 1, we can take k(l + x) + 1 w kx + 1 in Eq. (fT5|) . 
This case should be divided into further two cases, kx 2> 1 
or kx 1. We can take k(l + x) + 1 w fcx in the former case 
and fc(l + x) + 1 w 1 in the latter case. In the former case, 
we have 



\/eeh/e B (fc<l,fca;>l). 



(17) 



This is shown as case 4 in Table Q] In the latter case, we 
have 



07/1* 



,8/3 



(fc < l,fcx < 1). 



(18) 



This case should be divided into three sub-cases; (x <C 1), 
(z > 1 and kx 2 < 1), and (x > 1 and fcz 2 > 1). These 
are labeled as cases 1, 2, and 3, respectively in Table [T] We 
can calculate the cooling Lorentz factor 7 C for each case by 

Eq. ug. 

The ordering of the emission luminosities for various 
cases listed in Table [1] divides the GRB emission models 
into two groups. Since the synchrotron and SSC emission 
components can have very broad spectra, so that a con- 
dition Lph 3> max (L syn , L ssc ) is necessary for the photo- 
spheric emission to be dominant in the MeV energy range 
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Table X. Ordering of emission luminosities for various cases of Q = (t?/»7*) 8 / 3 , £ = € d e e h, and B = £d e B- 



Case Q 7c(or fc'),x Luminosities £ and B 



1 77 < V*(S < 1) k < 1, kx < 1 a; < 1 L ph > L up > L syn > L SS c 5 > max(£, B) 

2 x > 1, kx 2 < 1 L up > L ph > L ssc > L syn f > S > (£ 2 B) 1 /3 > B 

3 x > 1, kx 2 > 1 L up > L ssc > L ph > L syn £ > (£' 2 B) 1 /3 > g > 

4 k < 1, fcx > 1 L SS c > L up > L syn > L ph £ > (SB) 1 / 2 > g > B 

5 fc > 1 x > 1 L ssc > L syn > L up > L ph £• > B > Q 

6 x < 1, kx > 1, fcx 2 > 1 L syn > L ssc > L ph > L up B>£>£ 2 /B>g 

7 fcx 3> 1, kx 2 <C 1 L ayn > Lph > Lssc > L up B > £ > S > £ 2 /B 

8 a: < 1, kx < 1 L ph > L syn > L up > L ssc B > G > £ 



9 r; > r;«(g > 1) fc' < 1, k'x < 1 x <1 L ph > L up > L syn > L ssc G > max(£ - , B) 



rather than the synchrotron or SSC emission components. 
Thus cases 1, 2, and 8 are included in the photospheric emis- 
sion models, while cases 3, 4, 5, 6, and 7 are included in the 
synchrotron-SSC emission models. In this paper we focus on 
the former cases, and in particular cases 1 and 2, for which 
the UP emission component is dominant in the high-energy 
range. We will see that cases 1 and 2 can be consistent with 
the spectra of some LAT GRBs in Section [5] Here we show 
the description of the cooling Lorentz factor in cases 1 and 
2, 



7c ■ 



3m e TZ 



4m p ri,i(r;/7)*) 8 / 3 

2/3 



»» ,S (Tf) 



7/3^,2/3, ( a L \ 



(19) 



The UP luminosity is simply written as 



(20) 

Figure [2] shows an example of a broadband spectrum of 
the emission from the GRB jet for case 1. The parameters 
are taken as L 53 = 3,773 = 3, r a , 7 /T a = l,Hi = 2, /3 ph = 
—2.5, t v _2 = 1, P = 2.3, e d e e = 0.1, e d e B = 0.03, and 
a_L = 5. The source redshift is set to be z = 2. For these 
parameters, we have (?7/»7*) 8 ^ 3 — 0.4 > max^e^/i, ed£s), 
and r p h — 1 x 10 11 cm and rt ~ 2 x 10 14 cm, which satisfy 
our assumption 3> r p h (Eq.[8]). We generate the spectra of 
the four emission components using the approximate forms 
of the individual components. The photospheric emission is 
approximated as a smoothed broken power-law spectrum 



sF. 



A 



P h 



(e P h) 



1+ (eph) 



(21) 



where a constant A pn is given for the peak of eFf to be 
equal to L pa / (4ivd L ) (d,L is the luminosity distance) and 
we set s = 2. The UP, synchrotron, and SSC components 
are approximated similarly as smoothed broken power-law 
spectra with the eF e peak values being equal to L up /(47rd?,), 
is yn /(47rd| j ), and L ssc / {Aird 2 L ), respectively. The eF e spec- 
tral indices of the UP emission are given by 2 for e < £ up ,i, 
(3 - q)/2 for £ U p,i < s < e U p,h, (2 - p)/2 for s uPih < e < 
e up ,KN, and -p — 1 for e > £„ Pi kn, where £ up ,; ~ £ P h7i 2 > 
£ up ,h =^ £ P h7h (71 = min(7 m ,7 c ) and j h = max(7 m ,7 c )), 
£u P ,kn — (7/m e c 2 ) 2 /£ p h, and q = p (q — 2) for the slow 
cooling case 7 m < 7 C (for the fast cooling case 7 C < 7 m ). 
This is valid for the case of /3 p h $S —(2 + f>)/2, which is 
satisfied for our case /3 < —2.5 and 2 < p < 3. The spec- 



tr al indices of t he sy nchrotron and SSC emission are found 
in lSari fc Esinl (|200 ll ). Figure [2] shows that the overall spec- 
trum consists of a photospheric component in the mid-range, 
and synchrotron and UP components in the low and high 
energy ranges, respectively. We will use this approximate 
method to have a rough overall spectrum and to show that 
the emission from the photosphere and internal shock can be 
consistent with the time-resolved spectra of observed LAT 
GRBs, GRB 080916C, GRB 090510, and GRB 090902B in 
section 

The high-energy photospheric photons are absorbed 
by the e pair creation at the photosphere. Assuming 
that the photons are isotropic in the comoving frame 
of the jet, the opacity for the photons with the co- 
moving frame energy e' can be estimated by t^(e') ~ 
0.1<7T£L n w p h(r p h;£Lnn) r 'ph/(2?7), where £^ nn = m 2 c 4 /£', and 
"- P h( r ph;£') — iph(£'/ e P h)' 3ph /(47r'" p hCr7 2 £' p h) is the comov- 
ing photon density per unit energy around e' at the photo- 
sphere (cf., iLithwick fc Sarill200ll ). Here we only consider 
the high-energy part of the photospheric emission (with 
the index /3 pn ) as target photons, since the photon density 
~ e'n' ph (r p h; e') peaks at e' ph . If a calculated main target 
photon (comoving) energy for the pair creation break en- 
ergy is smaller than e' ph , the region is optically thin for the 
pair creation. The pair creation break energy, defined by 
T^(e') = 1, is given by 



-ph 



m e V 77* 



£ P h \" 2 -^ 



\ m e c 2 / 



(9xl0 3 )^(V) 8/ Vr/f ph+3 



(t 



£ph 



MeV 



-2-/3 ph 



GeV. 



(22) 



For the spectral model of the photospheric emission, we as- 
sume a spectral cutoff at this energy. 

The e pair creation opacity for the high-energy UP 
photons at the internal shock region can be estimated simi- 
larly. The main target photons for the high-energy UP pho- 
tons may be either the photospheric emission incident into 
the internal shock region or the UP emission itself. Since the 
incident photospheric emission is highly anisotropic at the 



8 K. Toma, X.-F. Wu, P. Meszdros 
10~ 5 f 1 ' 




e [MeV] 

Figure 2. Example of the model spectrum of case 1 (see Table [TJ. The spectrum consists of the photospheric emission (dashed line), 
the UP emission (dot- long-dashed line), and the synchrotron emission (dot-short-dashed line). The SSC emission is not shown since it 
is too dim. The parameters are L53 =3,773 = 3,r a ^/r a = = 2,/3 p h = — 2.5, t„,_2 = l,p = 2.3,£de e = 0.1, e^B = 0.03, = 5, 

and 2 = 2. The characteristic quantities are rj, ~ 4.3 X 10 3 , L ph ~ L(r)/r)* ) 8,/3 ~ 1 X 10 53 erg s _1 , e ph ~ 4fcT a (»7/»7*) 8/ ' 3 ~ 4 MeV, 
e5J 7 ~ 5 GeV, L up ~ Le d € e h ~ 2 x 10 52 erg s _1 , £ up , m ~ £ p h7m — 20 MeV, e up , c ~ e p h7? — 5 GeV, e up , K N — (r]m e c 2 ) 2 /e ph ~ 600 GeV, 
^syn — Ed € B ^up /fr/V) 8 / 3 ~ 1 X 10 51 erg s \ e S y n ,c — £syn(7c) — 80 eV, and e S yn,a — 60 eV. The intra-source pair creation opacity 
is estimated to be less than unity. We do not take into account the absorption by the pair creation through propagating in the 
intergalactic medium. The two-dot line describes a low-energy spectrum with an index a p ^ = —1.0, typical of observed GRBs, whose 
possible origins are discussed in section [6] 



internal shock region and the collision angle of the two pho- 
tons is typically very small, so that the pair creation opac- 
ity will be si gnificantly reduc ed from th at with the isotropic 
assumption ( Zou et all l201ll ) (see also iGranot et alj 120081 ; 
lAckermann et aL 2011). Yet, for simplicity, we estimate the 
pair creation opacity by the interaction with the photo- 
spheric emission with the isotropic assumption, which pro- 
vides a possible largest opacity. On the other hand, the 
target UP photons are much less anisotropic. We derive 
a minimum possible value of the break energy, which can 
be obtained as the lower one of the break energies caused 
by pair creations with the photospheric and the UP pho- 
tons calculated under the isotropic assumption. The pho- 
ton densities per unit energy in the relevant energy range 
are given by n ph (n;e') ~ L ph (e' / e' ph )^ h / (4irri cr] 2 s'l h ) 
and n' up (n;s') ~ L up (e' /^ p )" (1+9)/2 /(47rr 2 o)V 2 p , h ) for 
the photospheric and UP emission, respectively. Thus the 
possible minimum break energy is estimated by £™ 7 n ~ 
min(e™ ph , £™ up ) and 



"77 >ph 



40 



7c n 



1 2/3 ph +3 



\m e c z J 



-2-/3 ph 



x m e c 



1 x lor 

£ph 



(t 
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1 2-q 



7c XT] 



40 
2 x 10 1 

U GeV 



V m e c 2 J 
(4xl0- 3 ) 2 -(2£)" 



2-q 

XT) 3 



TeV, 



(23) 



where x is given by Eq. (|18[) . The opacity of the interac- 
tion with the photospheric emission (the UP emission) for 
all the high-energy photons is less than unity if £™ 7 " p h ^ 
(»7 2 m e c 2 ) 2 /£ ph (£™ up > (J7 2 m e c 2 ) 2 /£ up ,i). 

The energy £ syn ,a below which the synchrotron self- 
absorption effect is significant can be estimated by equat- 
ing synchrotron flux to the blackbody flux of the charac- 
teristic electrons in the shocked region, F E ^„, a = [(1 + 
Vdi]27rm e 7ch(£ sy n,a/ft) 2 (r 2 /rj) (e.g., iMeszaros fc Reesl 



1 1997^) . The characteristic Lorentz factor y c h of leptons is 
given by 7 a whose synchrotron energy is e S yn,a in the case of 
7a > Ji and otherwise by 7;. We show a formula of the ra- 
tio of £ s yn,a tO £ s yn,a = £syn ("fh), by which we can calculate 
£ S yn,a in a typical case "f a > Ji, 



£syn,a 



9 o T B 



iTi 



(24) 



where w = 2/(q + 4) for 7; < j a < jh and w — 2/(p + 5) for 
7a > "fh- The value of £ syn ,a is important for estimating the 
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observed synchrotron flux at the optical frequency (see the 
above example shown in Figure [2]). 

4.2 Case of 77 > 77* 

In this case the kinetic luminosity that can be dissipated in 
the internal shock is w L{rj/rj t )~ 1 . The internal shock ra- 
dius is estimated to be r, — 2ct v n 2 a~ L 2 , and our assumption 
Ti ^> rj reduces to 

<C 2ct v r/*a'£ 2 ~ 6 x 10 11 tv-i^^a^ cm. (25) 

1 a 

The minimum injection Lorentz factor of the leptons is 
given by Eq.©. The magnetic field energy density is U' B = 
LtdtB I '(47rr ^07777*), and then the characteristic synchrotron 
energy of leptons with Lorentz factor 7 is given by 



,(7) 



3heB' 2 
7 V* 



5 x 10" 1 7 2 L 



1/8 / r a ,r \ 

53 VrJ 



3/8 



X t~ 



-2 U L 



1/2 



eV. 



(26) 



V 0.1 

The photospheric, UP, synchrotron, and SSC luminosi- 
ties are estimated to be L p h ~ L, L up ~ xL, L syn ~ k'xL, 
and I/ssc — fc'a; L, where x is calculated by Eq. (ll2|) . and we 
have defined 



k' 



(27) 



(28) 



The cooling Lorentz factor of leptons is 

3m e 7Z 1 
7c ~ 4m p Tv (77/77*) fc'(l + a;) + f 

where the optical depth at r< can be written as Tj,j = 
( 7 "ph/ 7 "i)('7/ 7 ?*) -2/ ' 3 - The equations for the luminosities and 
7 C are the same as those for the case 77 < 77* by replacing 
(v/v*) 8 ^ 3 by 77/^7* and k by k! , so that the same argument 
for dividing cases for the order of the four luminosities can 
be made. In this case, however, 77 > 77* leads to k' < 1 and 
ed£eh/(r]/r)*) < 1. Then we only have the case of fc' <C 1 
and fe'a; <C 1, which we include as case 9 in Table[T] In this 
case the MeV and high-energy emissions may be dominated 
by the photospheric and UP components, respectively, sim- 
ilar to case 1. (The case of k' -C 1 and k'x 2> 1 leads to 
x ~ \J e e h/tB (2> 1), which is not consistent with the con- 
dition k'x 2 = e<ie e /i/ (77/77*) < 1.) The cooling Lorentz factor 
is calculated by 

3m e TZ ^ 1n 2 rl/4 / »V7 \ ~ 5/4 5/4 , -2 , on N 

7e " 47^(77/77*) " 7x10 (77) ^- 2ai - (29) 

The parameter a; and the UP luminosity are given by 
tdteh 



v* 



edteh. 



(30) 



(31) 



The pair creation break energies of the photospheric 
and the UP emission can be estimated by Eq. (|22|l with 
replacing (77/77* ) 8/3 by (77 / 77*)~ (8/3)(1+/3 p h) and Eq. |(23j with 
replacing 77 by 77* and setting x to be Eq. (|30[) . respectively. 
The synchrotron self-absorption energy can be estimated by 
the same equation as Eq. (I24|) . 



5 CASE STUDIES 

The previous sections have provided a general formulation 
of the emission from the photosphere and internal shock 
of the GRB jet. Here we focus on the cases in which the 
emission in the MeV energy range and in the high-energy 
range are dominated by the photospheric and the UP emis- 
sion components, respectively, instead of the synchrotron 
or SSC emission components (i.e., cases 1, 2, and 9 in Ta- 
ble [T]). We then show that such cases can be consistent with 
the observed time-resolved spectra of three of the best ob- 
served Fermi/LAT GRBs, especially in the energy range at 
and above the MeV spectral peaks. The detailed analysis re- 
sults of the observed spectra of brightest LAT GRBs, GRB 
080916 C, GRB 090902B a nd GRB 090510 have be en pub- 
lished (|Abdo et al.ll2009bl ld; lAckermann et al.ll2010aP l. and it 
is to these data that we apply our formulation. We assume 
that essentially the parameters W and ct v are in the efficient 
scattering regime (see section [3]), unless otherwise stated. 



5.1 GRB 080916C 

This burst is a long GRB that occurred at a redshift z ~ 4.35 
(corresponding to di ~ 1.2 x 10 29 cm). The spectral anal- 
ysis shows that all the spectr a of the five time- bins can 
be fitted by Band functions jAbdo et all l2009bl V These 
time-binned spectra are shown by the thin lines in Fig- 
ure [3] The high-energy spectral index of the second time- 
bin 3.6 — 7.7 s, /3 ~ —2.2, is significantly larger than that 
of the first time-bin 0.0 — 3.6 s, /3 ~ —2.6. This spectral 
hardening corresponds to the observed delay of the onset of 
the LAT emission with respect to that of the GBM emis- 
sion. The delay timescale in the cosmological rest frame is 
~ 5/(1 + z) s ~ 1 s. The high-energy spectral index seems 
stable after the second time-bin. We will see that the first 
time-bin spectrum is consistent with having only a photo- 
spheric component, while the subsequent time-bin spectra 
can be modeled as a photospheric plus a UP component, 
which mimic a smooth Band function. 

We first find the parameter values appropriate for the 
second time-bin spectrum. In our adiabatic jet model, we 
have a constraint /3 p i, < —2.5 (see Section [2}. Thus the 
observed hard high-energy spectrum of this time-bin with 
/3 ~ —2.2 has to be made by a combination of the pho- 
tospheric and UP emission. The spectrum can be fitted in 
this way as shown by the thick dashed line in Figure [3] The 
following conditions are required for producing such a spec- 
trum: L ph ~ 4.5 x 10 53 erg s _1 , s ph ~ 6.3 MeV, f3 ph ~ -2.5, 
L up ~ 1.1 x 10 53 erg s" 1 , e up , c ~ 2.5 GeV, e up , m < 42 MeV 
(which correspond to j c ~ 20 and 7 m < 2.6), and p ~ 2.8. 
Significant deviations from these conditions lead to a bumpy 
spectrum and/or deviate from the observed data with errors 
taken into account (see below for other conditions). The con- 
ditions indicate the slow-cooling case of the electron distri- 
bution, i.e., 7 m < 7 C . The fast-cooling case is not favored 
since it has a harder spectrum at e up , c < e < £ up , m which 
makes a dip in the spectrum. The ratio of the two lumi- 
nosities is Lu P /iph = x ~ 0.2, which indicates 77 < 77* for 
a reasonable value of £de e < 0.1 for the internal shock (see 
Equations 1181 and I30p . Thus this spectrum corresponds to 
case 1 (see Table [TJ. Equations for case 1, ©, (f20)) , ©, and 
(|19|l . translate the above conditions into five constraints on 
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Figure 3. Model fits of the time-resolved spectra of GRB 080916C. The thin lines are the observed data, and the thick lines are the 
model fits (solid lines are for 0.0 — 3.6 s, dashed lines are for 3.6 — 7.7 s, dot-long-dashed lines are for 7.7 — 15.8 s, dot-short-dashed lines 
are for 15.8 — 54.7 s, and the two-dot lines are for 54.7 — 100.8 s). Adopted parameter values are listed in Table[2] The spectral portions 
at and above the peak energies are well described by the model. The observed spectral parts below the peak energies are significantly 
softer than the model. This issue is discussed in Section [6] 



the model parameters, 




We can constrain the parameters L, r], r a /T a , TZ, and 
t v through these equations, by choosing reasonable values 
e^te ~ 0.1 and cll ~ 5 for internal shock. The parame- 
ters /3 p h and p have already been constrained directly by 
the model fit. The remaining parameter e^es does not af- 
fect the photospheric and UP spectra and it is constrained 
by the synchrotron emission contribution in the low en- 
ergy range (see below). The above first three equations (|32l 
El and [34ll provide L 53 ~ lift -1 , (r)/r)»f /3 ~ 0.41ft, and 
r a ,7/T a ~ 0.987i 3/2 . We have some allowed ranges of pa- 
rameters depending on the value of 7 m . For j m ~ 2.6, we 
have h ~ 0.98, and Equations (I35[) and ()36[) lead to H ~ 3, 
rj3 ~ 4.9, and t v ,-s ~ 0.2, where Equation ^ is used. If we 
adopt a smaller value, e.g., 7 m ~ 1, we have h ~ 0.45, lead- 
ing to 1Z ~ 8, ?73 ~ 7.4, and t Vt -2 ~ 0.02. The synchrotron 



emission fluxes, which are dominant below 10 keV, should 
be smaller than the observed fluxes around 10 keV, which 
provide upper bounds e^es < 0.1 for both cases. A larger h 
case has a smaller luminosity budget L. Comparisons of the 
prompt emission with the afterglow in GRBs generally imply 
high radiation efficiencies of the prompt emission, e 7 > 0.5 
(e.g.. |Panaitescu fc Kumari lioO^ . Thus we favor the case of 
the allowed smallest luminosity, i.e., the case of 7 m ~ 2.6, 
which we list the constrained parameter values in Tabled 

Next we consider the first time-bin spectrum. If the in- 
ternal shock for this component produces an electron en- 
ergy distribution with p ~ 2.8 similar to the second time- 
bin, the UP emission would have a hard photon index 
— (p + 2)/2 ~ —2.4 in the GeV energy range. Thus the 
UP emission has to be very dim compared with the pho- 
topheric emission, and the overall spectrum would consist 
mainly of the photospheric emission. The luminosity ratio 
should be I/ up /I/ p h = x < 0.06. The internal shock may have 
edte ~ 0.1 similar to the second time-bin, so that we require 
7/ > 77* (see Equations 1 181 and !30|) . Thus this spectrum corre- 
sponds to case 9 (see Table [T]). As shown by the thick solid 
line in Figured we can fit the spectrum mainly by the pho- 
tospheric emission whose luminosity, peak energy, and high- 
energy spectral index are set to be L p h ~ 3.1 x 10 53 erg s~\ 
e p h ~ 2.3 MeV, and /3 p h ~ —3.0, respectively. Equations 
(J3| determine the parameters L and r a /Y a by I/ P h and e p h, 
L53 ~ Lph,53 ~ 3 and r a ,7/T a ~ 20. The allowed ranges of 
r\, 1Z, and t v are broad since the UP spectrum is not tightly 
constrained. Equation (|3ip means that larger r\ decreases the 
UP luminosity. The model fit shown in Figure[3]is a case with 
the maximum contribution of the UP emission, i.e., a case of 
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Table 2. Best fit model parameters for GRB 080916C, GRB 090902B, and GRB 090510 



GRB 080916C^ 


Time-bin 


^53 


r a j/T a 


V3 




/3 p h 


t v ^ 2 


P 


e d e e 


e d e B 




v/v* 




0.0 - 3.6 s 


3 


20 


3.2 


2 


-3.0 


< 1.0 


2.8 


0.1 


< 1 




1.6 






3 


20 


3.1 


2 


-2.7 


0.1 


2.8 


0.01 


< 1 




1.6 




3.6 - 7.7 s 


10 


0.9 


1.9 


3 


-2.5 


0.2 


2.8 


0.1 


< 0.1 




0.71 




7.7 - 15.8 s 


5 


1 


3.2 


3 


-2.5 


0.6 


2.8 


0.1 


< 0.1 




0.61 




15.8 - 54.7 s 


3 


1 


2.5 


3 


-2.5 


0.7 


2.6 


0.1 


$ 0- 1 




0.52 




54.7- 100.8 s 


0.7 


2 


1.6 


3 


-2.5 


1.5 


2.6 


0.1 


$ 0.2 




0.55 


GRB 090902B a 


Time-bin 


^53 


r a ,r/T a 


13 




/3 ph 


t v ,-2 


V 




£d£B 




v/v* 




4.6 - 9.6 s 


16 


2 


3.5 


3 


-4.0 


10 


2.9 


0.1 


0.1 




0.52 


GRB 090510" 


Time-bin 


£53 


r a ,l/Y a 




Til 


/3 ph 


t v ,-2 


P 




£d£B 




v/v* 




0.5 - 0.6 s 


1 


3 


73 


2 


-4.8 


0.3 


2.3 


0.1 


< 1 




26 






1 


3 


5.4 


2 


-4.8 


0.1 


2.3 


5 x 10~ 3 


< 1 




1.9 




0.6 - 0.8 s 


1 


0.3 


4.7 


2 


-3.0 


0.2 


2.3 


0.1 


< 7 x 10" 


-3 


0.72 




0.8 - 0.9 s 


20 


2 x 10" 3 


3.8 


0.5 


-2.8 


0.2 


2.05 


0.15 


< 3 x 10" 


-6 


0.18 






10 


0.01 


1.6 


0.6 


-2.5 


7 


2.15 


0.15 


< 4 x 10" 


-6 


0.12 



a These values are example sets of parameter values for producing the fitting spectra, for which L is the allowed smallest values 
for each time-bin. See texts for detailed explanations how to specify the other values. Two sets of the parameter values are shown 
for the time-bins of 0.0 - 3.6 s of GRB 080916C, 0.5 - 0.6 s and 0.8 - 0.9 s of GRB 090510. The parameter a L is taken as 5 for all 
the time-bins except that we take = 2 for the second set for the 0.0 — 3.6 s time-bin of GRB 080916C and for the second set for 
the 0.5 - 0.6 s time-bin of GRB 090510. 



the minimum 77 for the fixed values e^ee ~ 0.1 and p ~ 2.8. 
This model fit requires the luminosity and characteristic en- 
ergies of the UP emission to be L up ~ 2.0 x 10 52 erg s _1 , 
£u P ,h ~ 27 MeV, and e uPl ; < 10 MeV, respectively (corre- 
sponding to 7/j ~ 3.4 and 7; < 2.1). Equations ([9]) and (|29|) 
indicate that we can have similar values of the model pa- 
rameters 1Z and t v as for the second time-bin, i.e., we have 
jm ~ 3.4 and j c < 2.1 for TZi ~ 2.4 and t„ _ 2 < 1.0. Since 
7c < 7m, we have h = 1, so that Equation ()31[) provides 
773 ~ 3.2. 

Alternatively, we may consider a case of ede e <C 0.1 for 
the first time-bin. For the case of 77 > 77*, the final Lorentz 
factors of the shells are given by 77* , which may not be much 
different, i.e., a_L < 2, since it depends weakly on L, r a /Y a , 
and 1Z. Then the internal shocks may only cause a weak 
dissipation, with e^e <C 0.1. The observed first time-bin 
spectrum can be fitted by a photospheric emission only, with 
/3 p h — —2.7 in this case. The model parameters L and r a /Y a 
are constrained to be the same as the above case e^e ~ 0.1. 
For edte ~ 0.01 as an example, we can find a set of parameter 
values 77, 1Z, and t v similar to those for the second time-bin, 
which are shown in Table [2] 

The spectral shapes of the other three time-bins are 
very similar to that of the second time-bin. They can also 
be produced by a combination of the photospheric and UP 
components in our model, where the luminosity ratios of 
the two components are required to be relatively large, i.e., 
x ~ 0.3, which indicate case 1, similar to the second time- 
bin. We can constrain the model parameter values for the 
three time-bin spectra in a similar way to the second time- 
bin spectrum. The parameters constrained for the case of 
the allowed smallest luminosity budgets are summarized in 
Table [2] The rightmost column represents the value of 77/77* 
calculated for the adopted parameter values. 



The e pair-creation break energies under the isotropic 
photon field assumption is given by eJL m i n /(l+ z ) ~ 10 GeV 
for the first time-bin and e~~ m i n /(l + z) > 50 GeV for 
the other time-bins. Since the target photon field is highly 
anisotropic, the real break energies may be much larger than 
those values (see references listed above Equation I23[) . so 
that we have neglected the pair-creation breaks in Figure 
We confirm that the second-order UP emission is negligible 
for each time-bin. For the second to fifth time-bin spectra, 
we have 777 c m e c 2 7 2 e up , c , so that the KN effect suppresses 
the second-order UP emission. For the first time-bin, the 
second-order UP emission flux is just much smaller than 
the other emission components. The parameters satisfy our 
assumptions of r, 3> 77 for the first time-bin (Eg. I25[) and 
t% 3> 7-ph for the other time-bins (Eq. [8]). 

We have shown that the parameter regime should shift 
from 77 > 77, into 77 < 77, to reproduce the observed first and 
second time-bin spectra in our model. This shift is related 
to the large decrease of r a /F a from ~ 2 x 10 8 cm in the first 
time-bin to ~ 1 x 10 7 cm in the second time-bin. At this 
transition, the luminosity ratio of the photospheric and UP 
components x — L up /L p h (Equations 1301 and I18[) increases 
as TafYa decreases (and then 77* increases) and possibly tdte. 
increases because of the increasing differences of the final 
Lorentz factors of the shells. As we discuss in Section [6] 
the simulations of th e jet dynamics in the progenitor star 
l|Morsonv et al.l l2007h suggest that r a /T a may undergo a 
sudden decrease as time progresses. Thus this transition may 
be interpreted as the reason for the observed delay of the 
LAT emission onset. After the transition, the parameter sets 
constrained for the allowed smallest luminosities indicate 
that L, 77, and t v evolve monotonically, while r a /Y a , 1Z, /3 p h, 
and p are stable (note that the values of eae e and a_L are 
given to determine the other parameters). 
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Figure 4. Model spectra for the second time-bin spectrum of GRB 080916C, for different values of t v and 1Z and fixed values of L, 
Ta/^a, (t?/'?*) 8 / 3 ! /3 p hi P> e d e e an d e d e s as shown in Table[2] Left: The thick solid line is the model spectrum using the best fit parameter 
values shown in Table [2] except for taking t' v = 2t v , which is comp ared with the thin lines of the Band function fitting the observed data, 
with the dashed lines showing the la errors jAbdo et al.ll2009bri . Right: The thick solid line is the model spectrum using the best fit 
parameter values shown in Table [2] except for taking TV = 1Z/5, which is compared with the same observed data as the Left panel. In 
both panels the thick dashed lines and the thick dot-dashed lines represent the photospheric and UP components, respectively. 



The variability timescale is as small as t v ~ 10 -3 s, 
which is larger than the light crossing timescale across a 
Schwarzschild black hole of mass ~ 10 M e , ~ 1(T 4 s. Such 
a timescale (t v (l+z) ~ 5 ms in the observer-frame) may not 
be resolved by the 7-ray detectors that observed this burst 
as far as we know. The most sensitive detector among them 
for the 100 keV- 30 MeV range is INTEGRAL, which has a 
resolution down to 50 ms and found a variability timescale 
~ 100 ms (IGreiner et al. 2009). An unsteady mass accretion 
from a convectively unstable torus onto the central com- 
pact ob ject might produce such a large timescale of the flux 
change (|Sekiguchi fc Shibatdl20ich . or the interaction of the 
jet with the stellar envelope could induce a variability with 
light crossing timescale from the side to the axis of the jet 
~ R,6j(l + z)/c ~ 100 (R,/5x 10 10 cm)(6» J /0.01) ms, where 
R* is the stellar radiu s and 9j is the opening angle of the jet 
|Morsonv et alj|2010h . 

As seen in Figure(3] the low-energy portions of the spec- 
tra of the photospheric emission are much harder than the 
observed spectral portions below the peak energies. This 
is currently a generic problem for all photospheric mission 
models, which will be discussed in Section [5] 

Interestingly, the observed peak energies and peak eF e 
fluxes in the second to fifth time-bins are clearly consistent 
with the famous pe ak energy-luminosity relation (so-called 
Yone t oku relation; lYonetoku et al.l |2004| ; iGhirlanda et al.l 
120051 . 120101 ). which implies e p h oc L^ 2 in our model. This 
translates into a relation r\ oc ^/"Wr^/^^byEqua- 
tions H}. Considering that r a /T a and 1Z are not much 
changed in our results and the dependences in the relation 
are weak, the dominant relation should be 7/ oc L 7 / 16 (or 
roughly 77 oc L 1 ^ 2 ). A degree of fine tuning appears needed 
to satisfy this relation; this is similar to other GRB emission 
models, whic h also require fine tunings of the parameters of 
the jet (e.g., IGhirlanda et all feoid : iRees fc Meszarosll2005l : 
iThompson et al.ll2007l : llokal^OlOT ) . 

Furthermore, we should clarify how tightly the param- 
eter values related to the UP component are constrained 
for reproducing the smooth Band-like spectrum by a su- 



perposition of the two components for each time-bin. Tak- 
ing the second time-bin as an example, we make differ- 
ent model spectra by fixing L, r a /F a , (r//r] t ) 8 ^ 3 , and /3 p h 
to have the same photospheric component, while changing 
the other parameters to see whether we obtain a smooth 
Band-like spectrum. The left model in Figure [4] is the re- 
sult for a value of t v that is 2 times larger than the value 
in Table [2] which is marginally consistent with the data 
at la level. Even larger t v , which leads to larger j c (and 
then larger £ up , c and smaller h), makes the UP component 
dimmer and then the overall spectrum significantly devi- 
ates from the data. We also examine models for smaller 1Z, 
which leads to smaller f c (with larger 7m). This case keeps 
h ~ 1 (and then L up = Ltdteh ~ const.), making the overall 
model spectrum somewhat bumpy. For a value of TZ that is 
2 times smaller than the value in Table [2] the bumpy model 
spectrum is marginally consistent with the data at la level 
at e/(l + z) < 1 GeV but does not reproduce the flux at 
e/(l + z) > 1 GeV because of smaller £ up , c . Even smaller 
1Z leads to the fast-cooling case j c < y m . The spectrum be- 
low the UP peak energy is harder than the slow-cooling case 
with p ~ 2.8, so that the model spectrum is more bumpy, 
although a higher £ up , m can explain the emission at the high- 
energy range. The right model in Figure [4] is the result for 
a value of 7Z that is 5 times smaller than the value in Ta- 
ble [21 which is still marginally consistent with the data at 
la level. For a smaller 1Z, however, the overall model spec- 
trum significantly deviates from the data in the middle range 
~ 10 MeV — 1 GeV. This parametric study indicates that 
the parameter spaces of t v and 1Z are limited within a fac- 
tor of ~ a few for obtaining a Band-like function by the 
superposition of the two components. 

This means that in order to obtain the Band-like spec- 
tra through the second to fifth time-bins, we need significant 
fine tuning of the model para meters t v and 7Z (as w ell as /? p h 
and p). However, we note that lAbdo et all |2009bh conclude 
"Compared to the null hypothesis that the data originated 
from a simple Band GRB function, adding the additional 
power-law component resulted in a probability of 1% that 
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Figure 5. Model fit of the spectrum of GRB 090902B in the 4.6 — 9.6 s time bin. The thin line is the observed data, and the thick 
line is the model fit. Adopted parameters are listed in Table [2] The photospheric and UP emission components are consistent with the 
observed spectra of the Band component at and above the peak energy and the distinct high-energy power-law component, respectively. 
The observed low-energy spectrum of the Band component is softer than the model. This issue is discussed in Section [6] A deviation 
from the Band component at < 40 keV could be explained by a contribution from the synchrotron component. 



there was no additional spectral component for this (fourth) 
time bin." Taken at face value, this would imply that the 
fourth time-bin data may have an additional high-energy 
spectral component besides the Band function although the 
significance level is l ow (~ 2a level). Fur thermore, the light 
curve at < 10 s in lAbdo et all J2009bl ) shows the abrupt 
steepening break of the GBM light curve at the fifth time- 
bin while the LAT flux decays st ably, which implies a two 
component origin at this time-bin. IZhang et al.l (|2010f ) claim 
that the background uncertainty significantly affect the re- 
sults of the spectral analysis due to the low count rate at 
such late times. These implications increase the plausibility 
of our two component model. 



5.2 GRB 090902B 

This burst is a long GRB that occurred at a redshift z ~ 1.82 
(which corresponds to dh — 4.3 x 10 28 cm). The r e sults o f 
the time-resolved spectral analysis by lAbdo et all (|2009ci ) 
show that all the spectra can be fitted by a Band function 
plus a distinct power-law function. They did not list the 
flux normalizations of the two components and only show an 
overall spectrum for the second time-bin 4.6—9.6 s, which we 
show by the thin line in Figure [5] We apply our model only 
for this time-bin spectrum. The LAT emission starts ~ 3 s 
after the GBM trigger within the first time-bin 0.0 — 4.6 s. 
The delay time in the cosmological rest frame is ~ 3/(1 + 
z) s ~ 1 s, which is similar to that of GRB 080916C. 

The Band and distinct high-energy power-law compo- 
nents of the second time-bin spectrum can be straightfor- 
wardly modeled by the photospheric and UP emission, re- 



spectively, as shown by the thick line in Figure [S] The ob- 
served data of Fermi/GUM show a clear deviation from 
the Band function at < 40 keV, and the thin line is ob- 
tained by assuming that the low-energy excess and the high- 
energy emission a t > 50 MeV are the same power-law emis- 
sion component (|Abdo e t al. 20093). I n °ur model, how- 
ever, the low-energy excess could be a contribution from the 
synchrotron component, which is separate from the high- 
energy UP component. If the superposition of the multi- 
ple photospheric emission from different shells can repro- 
duce the observed low-energy power-law portion of the Band 
component (see Section [5] for more discussion) and it ex- 
tends even below ~ 40 keV, the contribution from the syn- 
chrotron component could explain the observed low-energy 
excess. The model fit is obtained by setting the conditions 
L ph ~ 2.8 x 10 53 erg u" 1 , e ph ~ 2.4 MeV, p ph ~ -4.0, 
L up > 1.3 x 10 52 erg s -1 , e up , c > 40 GeV, e up , m < 16 MeV, 
and p ~ 2.9 (which correspond to 7 C > 130 and 7 m < 
2.6). For these values we have L up /L p h = x > 0.06, and 
h = (7c/7m) 2_p /(3 -p) < 0.2. These lead to r? < 77, for 
a reasonable assumption of e^e < 0.1, and thus this spec- 
trum corresponds to case 1, similar to the second time-bin 
of GRB 080916C. From these conditions and given values 
e^ee ~ 0.1 and ~ 5, we can constrain the model param- 
eter values, in a similar way to that for GRB 080916C. A 
larger h (for a smaller 7 C and a larger r j m ) leads to the 
smaller luminosity L, while in order to explain the low- 
energy excess as the synchrotron emission at e > £ sy n,c with 
a reasonable range e^es <■ 0.1, we require 7 C > 350 (for 
which L up > 1.6 x 10 52 erg s _1 is needed to explain the 
high-energy emission). In Table [2] we show the parameter 
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Figure 6. Model fits of the time-resolved spectra of GRB 090510. The thin lines are the observed data, and the thick lines are the 
model fits (the solid lines are for 0.5 — 0.6 s, and the dashed lines are for 0.6 — 0.8 s). Model parameter values are listed in Table[2] The 
observed spectral parts below the peak energies are significantly softer than the model. This issue is discussed in Section [6] 



sets constrained for the case of "/ c ~ 350, 7 m ~ 2.6, and 
Lu P ~ 1.6 x 10 52 erg s _1 . They satisfy our assumption of 
Ti S> r p h (Eq. [8]). The internal shock region is found to be 
optically thin for e pair creation. The second-order UP 
emission is confirmed to be negligible. 

It has been reported that the Band component of this 
burst can be fitted by a multi-temperature blackbody spec- 
trum <|Rvde et alj|2010h . The observed low-energy spectral 
slope of the Band component a ~ 0.07 is relatively hard and 
near (but still softer than) the single-temperature thermal 
model Q p h = 1, and the high-energy spectral slope /3 ~ —3.9 
is very steep. These strongly favor a photospheric origin of 
the Band component of this burst. 

The first time-bin of GRB 090902B includes a very 
small number of the high-energy photons, which implies 
larger 77/17*. Thus the LAT onset delay in this burst may be 
consistent with the interpretation that the parameter regime 
shifts from 77 > 77, in the first time-bin into 77 < 77, in the 
second time-bin, similar to the case of GRB 080916C. 

5.3 GRB 090510 

This burst is a short GRB that occurred at a redshift 
z ~ 0.90 (corresponding to d,L — 1.8 x 10 28 cm). This has 
a precursor ~ 0.5 s earlier than the main burst. The anal- 
ysis results of four time-binne d spectr a of the main burst 
are given in lAckermann et all (^OlOal ). The first and sec- 
ond time-bin spectra are shown by the thin lines in Fig- 
ure [6] and the third time-bin spectrum in Figure [7] The 
first time-bin spectrum can be fitted by a Band function 
only, while the second and third time-bin spectra can be 
fitted by Band plus distinct power-law functions. For the 
last time-bin 0.9 — 1.0 s, the emission was detected only in 



the LAT energy range (100 MeV — 2 GeV), whose spectrum 
is fitted by a power-law function. We will not fit this last 
time-bin spectrum with our model but we discuss a possible 
explanation for it in the context of our model. The onset 
delay of the distinct high-energy component with respect to 
the Band component in the cosmological rest frame is es- 
timated to be ~ 0.1/(1 + z) s ~ 0.05 s. It is remarkable 
that the distinct high-energy component is brighter than 
the Band component in the third time-bin, which is unique 
among all the observed LAT GRBs. 

We first consider the second time-bin spectrum. This 
may be explained by a combination of the photospheric and 
UP emission components. The model parameters for this 
spectrum can be found in a similar way to the case of GRB 
090902B. The luminosity ratio is required to be -I/ up /L p h = 
x > 0.1, and we estimate that h < 0.7. These lead to 77 < 77* 
for a reasonable assumption e^te < 0.1, and thus this spec- 
trum corresponds to case 1, similar to the second time-bins 
of GRB 080916C and GRB 090902B. The best fit model 
for this spectrum is shown by the thick dashed line in Fig- 
ure H for which L ph ~ 1.7 x 10 53 erg s" 1 , e ph ~ 9.7 MeV, 
/3 ph ^ -3.0, L up > 2.3 x 10 52 erg s~\ e up , c > 20 GeV, 
£up,m iS 44 MeV, and p ~ 2.3 are required. From these 
conditions and given values edt e ~ 0.1 and cll ~ 5, we can 
constrain the model parameter values in a similar way to the 
cases for GRB 0809 16C and GRB 090902B. The parameter 
sets constrained for the allowed smallest luminosity (i.e., for 
£u P , c ~ 20 GeV, e up , m ~ 44 MeV, which lead to the largest 
h, and L up ~ 2.3 x 10 52 erg s -1 ) are shown in Table(2] 

Next we consider the first time-bin spectrum. This 
should be produced by the photospheric emission only and 
the luminosity ratio is required to be L up /L p h = x < 
3 x 10~ 3 . For a reasonable value £d£ e ~ 0.1, we have 77 > 77,. 
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Figure 7. Model fits of the third time-bin (0.8 — 0.9 s) spectrum of GRB 090510. The thin line is the Band plus power-law function 
fitting to the observed data and the thin dashed lines represent the lcr errors. The thick solid and dashed lines are our model fits. Model 
parameter values are listed in Table [2] See texts for details. 



This spectrum corresponds to case 9 (see Table [T]). The 
model fit shown by the thick solid line in Figure [6] requires 
that the photospheric emission has L p h ~ 1.4 x 10 53 erg s -1 , 
e P h ~ 5.5 MeV, and /3 p h ~ —4.8. Equations © determine 
L53 ~ L p h,S3 ~ 1 and r a ^/Y a ~ 3. These lead to a very 
large 77/77* ~ Led£eh/L up ~ 30(e<je e /0.1), indicating a very 
large 77 ~ 10 . The parameters for the model fit shown in 
Figure [6] are shown in Table [2] 

We can also consider a case of 77 > 77* and edt e < 0.1 
for the first time-bin, similar to the first time-bin of GRB 
080916C. We show a result for a case of ede e ~ 5 x 10~ 3 
in Table [2l in which 77/77* ~ 2 and 77 ~ 5.4 x 10 3 . In this 
case Equation provides ■y m ~ 0.1 < 1, which means in 
reality that most of the leptons have 7 ~ 1 and only a small 
fraction ~ 0.1 of them participate in the non-thermal power- 
law acceleration. The UP component spectrum has a bump 
around £ pn , which makes a small contribution to the overall 
spectrum. 

Next we find the parameter values appropriate for the 
third time-bin spectrum. In this time-bin the UP luminos- 
ity has to be higher than the photospheric luminosity. This 
may correspond to case 2 (see Table [T]), unlike the cases 
considered above. The spectrum can be fitted as shown by 
the thick dashed line in Figure [Jj where the emission com- 
ponent at e/(l + 2) < 0.1 MeV is the SSC emission (the 
second-order SSC emission is hidden by the UP emission). 
The characteristic quantities required for this model fit are 
L ph ~ 1.8 x 10 52 erg s -1 , e ph ~ 3.6 MeV, /3 ph ~ -2.8, 
L up > 2.3 x 10 53 erg s~\ e up ,c > 100 GeV, e up , m < 50 MeV, 
and p ~ 2.05. For a somewhat larger edt e ~ 0.15, we can 
constrain the model parameter values in a similar way to 
case 1, and the results are shown in Table [2] The value 
7"a,7/r a ~2x 10 -3 is much smaller than that in the second 



time-bin. This leads to r a > 2 x 10 4 cm, which is too small, 
compared to a Schwarzschild radius ~ 3 x 10 5 cm of the 
central compact object of mass ~ Mq. If we take a smaller 
Lup and £ up , c which are consistent with the data at la level, 
i.e., L up ~ 1.2 x 10 53 erg s _1 and e up ,c ~ 25 GeV, we have 
r a /T a ~ 6 x 10 4 cm (and we have L53 ~ 10 smaller than 
the larger L up model). At such a late phase of the prompt 
emission, the external sh ock emission might con tribute to 
the high-energy emission (|De Pasauale et alJbOltfh but this 
poss ibility is suggested to be unlikely ( He et al. 2010T ) (see 
also lLiu fc Wariell201Ch . Below we show another sets of pa- 
rameter values which could explain the third time-bin emis- 
sion at la level together with the fourth time-bin emission, 
and have an even larger r a /Y a . 

Finally we discuss the fourth time-bin spectrum. In this 
time-bin, only the high-energy photons are detected, which 
requires much larger e^ee and/or much smaller r a . These do 
not seem realistic. Here we show an example of another pa- 
rameter set for the third time-bin spectrum, for which the 
duration of the UP emission can be larger than the pulse 
width of the MeV photospheric emission and stays bright 
even in the fourth time-bin without the MeV emission. Here 
we consider the inefficient scattering regime for the third 
time-bin (see Section [3] we have assumed the efficient scat- 
tering regime for the first and second time-bins, say a case 
of W ~ ctv), in which the widths of the two shells W are 
not much different from those for the earlier time-bins W, 
but the separation of the two shells ct v is much larger than 
those for the earlier time- bins ct v . In this regime we consider 
a case of aj^ct-u ~ W + ct v ~ 2cr5„, for which the leptons 
in the internal shock of the given two shells with W and t v 
up-scatter the photospheric emission from a third shell at a 
distance ct v behind these two shells. In this case the angular 
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spreading timescale of the UP emission is estimated to be 
~ 3a^ 1 t t ,, which is much larger than W/c and t v , so that 
we have 8t up ~ 3a^ t„ (see Eq.[7J|. Since 6t p h ~ W/c, which 
means that the UP emission can last much longer than the 
MeV emission. If 8t ap (l + z) is comparable or larger than the 
duration of the time-bin, 0.1 s, this scenario could explain 
the fourth time-bin spectrum. Below we show that the model 
fit of the third time-bin spectrum in this scenario provides 
a t v consistent with the above temporal conditions. 

We interpret the observed Band and distinct high- 
energy components in the third time-bin as the photospheric 
emission of the third shell and the UP emission produced 
by up-scattering of the photospheric emission of the third 
shell by the internal shock of the first and second shells, 
respectively. The model fit is shown by the thick solid 
line in Figure which is marginally consistent with the 
data at la level. The photospheric emission of the third 
shell has L ph ~ 1.8 x 10 52 erg s~\ e ph ~ 5.5 MeV, 
and /3 p h — —2.5. (Hereafter a tilde denotes a quantity 
of the third shell.) The UP emission is required to have 
L up ~ Le d e e h ~ 1.2 x 10 53 erg s" 1 , e up , c ~ e ph -y^ ~ 25 GeV, 
e uP:m ~ e p h7m ~ 180 MeV, and p ~ 2.15. The condition 
for L up gives us L53 ~ 10 for a given ed£ e ~ 0.15. The pho- 
tospheric flux of the second shell should be lower than the 
observed flux level. We take L p i r ~ 3.5 x 10 51 erg s _1 and 
e P h ~ 0.5 MeV which produces the spectral bump below 
the peak energy of the main photospheric emission in Fig- 
ure [7] and requires r a ,7 /r„ ~ 0.01 and (v/v*) S/3 ~ 4 x 10" 3 
through Equation Q. Equation (J9]) leads to IZi ~ 0.6. 
Then we have 77* ,3 ~ 13 and 773 ~ 1.6. Finally the re- 
maining parameter t v is determined by an equation for 
7c which is given by Equation (|19[) multiplied by a fac- 
tor Lph/iph ~ 0.2. This factor leads to t Vl —2 ~ 7, much 
larger than in the above simple two-shell scenario. The ob- 
served duration of the UP emission pulse is then estimated 
to be <5t up (l + z) ~ 3a£ i„(l + z) ~ 0.08 s. This is com- 
parable to the durations of the third and fourth time-bins 
and roughly consistent with the condition a~/ 2 ct v ~ 2ct v , so 
that the high-energy emission without a corresponding MeV 
emission in the fourth time-bin could be explained by the 
UP emission with large angular spreading time. 

The value r a /T a ~ 10 5 cm is comparable but still some- 
what smaller than a Schwartzschild radius of mass ~ Mq . A 
larger r a /Y a would lead to a brighter spectral bump of the 
photospheric emission from the second shell at ~ 0.1 MeV 
in Figure [3 as L ph oc (r a /T a ) 2/3 and e ph oc (r a /r a ) 1/6 in 
Eq. (|4]), violating the observed data. This might suggest that 
a large fraction of the jet ener gy at r a is not therm al but 
Poynting flux for this burst (cf. lzhang fc Pe'edl200gh . 

To summarize this scenario, the pulse widths of the UP 
component St up ~ W/c + 3a^t v are similar to those of the 
photospheric component 5i p h ~ W/c (as ~ 1 ms) at the first 
and second time-bins, while at the third time-bin the UP 
pulse width only increases significantly to St up ~ Sa^t v ~ 
40 ms, which may stay bright even at the fourth time- 
bin. This may be consistent wit h the results of the cross - 
correlation function analysis by lAckermann et all (|2010al ) 
that there is no correlated variability between the GBM and 
LAT emission(3 In this scenario there should be a large in- 



4 Nevertheless, we can see some narrow spikes of the LAT emis- 



terval t v (l + z) ~ 0.1 s of the photospheric emission before 
the third time-bin, which seems consistent with the observed 
quiescent time in the 260 keV — 5 MeV light curve. Such a 
variability timescale ~ 0.1 s might arise from the mass accre- 
tion o n the central c ompact object from an inhomogeneous 
torus (Rosswoe 20071) o r from the interaction of the jet with 
the dense environment (|Morsonv et alj|2010h . 

The parameter sets we found for the three time-bin 
spectra are summarized in Table [5] We have shown that the 
parameter regime should shift from 77 > 77, in the first time- 
bin into rj < 77, in the second time-bin, which is related to 
the large decrease of r a /T a and corresponds to the increase 
of the luminosity ratio x = L up /L p h and the delayed onset 
of the distinct high-energy component, similar to the cases 

of GRB 080916C and GRB 090902B. The second and third 

1/2 

time-bin spectra satisfy the relation e pn oc L ( , similar to 
the case of GRB 0809 16C, which requires a fine tuning of 
the parameters L, r a /T a , rj, and 7Z. 

The synchrotron and SSC fluxes should be lower than 
the observed flux levels, which put constraints on edtB, 
which is much tighter than the cases of GRB 080916C and 
GRB 090902B. In particular for the third time-bin, we need 
£de_B < 3 x 10~ 6 . A typical value ~ 0.25 leads to a 
constraint es < 10~ 5 . This is not implausible since some 
external shocks driven by GR Bs have been suggested to 
have such small values of en ([Panaitescu fc Kumar] |2002| ; 
iGranot fc Guettall2003l ; IPe'er fc Waxmanll2004r i. The inter- 
nal shock region is estimated to be optically thin for the 
e pair creation for the second and third time-bins, and we 
find £™ ph /(l + z) ~ 100 GeV for the first time-bin. The 
parameter sets for all the time-bins satisfy our assumptions 
of Equations (JSj and (125J1 . The second-order UP emission is 
negligible for all the time-bins. 



6 SUMMARY AND DISCUSSION 

GRB jets are thought to consist of many successive shells 
moving at relativistic speeds. These naturally lead to vari- 
able photospheric emission around the MeV energy range, 
and also to internal shocks above the photosphere. Energy 
dissipation near the photosphere caused by e.g., internal 
shocks, excited plasma waves, and/or interaction of the jet 
with the dense environment, or nuclear collisions between 
protons and neutrons is expected to make the photospheric 
emission having the power-law tail spectrum above the peak 
energy e pn , with the photon index /3 pn < — 2.5 under the as- 
sumption of the roughly adiabatic evolution of the jet (see 
Section [2}. We have generically studied the temporal and 
spectral properties of the radiation from the photosphere 
and the internal shock in the jet which is accelerated by 

sion correlated with GBM pulses throughout the whole prompt 
phase. This may not be inconsistent with our second (large 8t up ) 
scenario for the third time-bin, since the third time-bin model 
spectrum indicates that the high-energy tail part of the photo- 
spheric emission which has sharp variability contributes to ~ 50% 
of (~ 10% of) the observed photon numbers at the entire LAT 
range > 20 MeV (at > 100 MeV). Also, the UP emission comes 
from a large area of the shell with solid angle of < 3/T (see Sec- 
tion [3} , which could have angular inhomogeneity leading to the 
sharp variability at the both third and fourth time-bins. 
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the thermal pressure. We have not considered acceleration 
by processes related to magnetic fields. We have shown that 
the photospheric emission is efficiently up-scattered to the 
high-energy range by the electrons (and positrons) acceler- 
ated in the internal shocks, for the efficient scattering regime 
(V 2 3 /Vl)t v < W/2 which includes the typical case W ~ ct v 
(see Section [3]), when the radiation from the internal shocks 
consist of the UP, synchrotron, and SSC emission compo- 
nents. 

Our generic arguments show that the ordering of the lu- 
minosities of the emission components depends on the val- 
ues of (f?/'?*) 8 '' 3 , Edfe^, and £des (see Section [4] and Ta- 
ble [TJl. A condition (ri/r] t ) 8 ^ 3 2> max(ede e /i, e^es) is re- 
quired in order to have the photospheric and UP compo- 
nents dominant in the MeV energy range and in the high- 
energy range, respectively, rather than the synchrotron and 
SSC components, and L p h 3> L up as typically observed in 
LAT GRBs. (We can also have a case of L up 3> £ P h for 
a condition e d e d h > (v/v*) 8/3 > [{^e) 2 ^d^B] 1/3 , which 
has been applied to the third time-bin spectrum of GRB 
090510.) For reasonable values of e^e and e^es both < 0.1 
for internal shocks, that condition reduces to 77 > 77, ~ 
2.8 x W 3 Ll / 3 ' i (r a , 7 /r a y 1/4 TZ 1 1 /4 (Equation [2]). In this case 
the photospheric emission has luminosities L p h ~ L, and 
its peak energy is given by e p h ~ 8 L^ 4 {r a ,7 /T a ) _1/ ' 2 MeV 
(Equations [T] O and 2]), which can be consistent with the 
Band components of the observed LAT GRBs as well as 
other GRBs since we can have 10 _1 < r a ,7/r a < 10 2 in 
principle (but there are several issues to be addressed on 
the spectral shapes and their temporal evolution; see be- 
low). This implies that the GRB jets should typically have 
a relatively large bulk Lorentz factor at the prompt emission 
site, ~ 77*, in our model. 

The UP emission component is a good candidate for 
the observed high-energy emission of LAT GRBs in the 
prompt phase. The long-lived high-energy emission after the 
prompt phase may be re l ated to the external shock (e.g. . 



De Pasquale et al. 2010]; Kumar fc Barniol Duran 



Ghisellini et al.ll20 id : iHe et al Il2010l : iLiu fc Wandl20101 ). We 
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have performed case studies of the prompt emission of GRB 
080916C, GRB 090902B, and GRB 090510, and fitted their 
time-binned spectra by our model typically with reason- 
able sets of parameter values (see Section [S] and Tabled- 
For GRB 090902B and GRB 090510, the UP emission cor- 
responds to the observed distinct high-energy component, 
while for GRB 080916C, the UP emission is dominant in the 
high-energy range but the photospheric plus UP emission 
mimic the Band function between ~ 10 keV and ~ 10 GeV. 

We note that fine tuning of the parameter values ap- 
pears to be required for some spectra in addition to keep- 
ing V }t V*- in order for the combined photospheric and 
UP components to mimic a simple Band function fit of 
the spectra of GRB 080916C through the second to fifth 
time-bins, significant fine tuning of several parameters is 
required. Such a composite model is, however, compatible 
with the presence of an additional high-energy component, 
which is weakly suggested in the data of the fourth time- 
bin, and also the difference of the temporal behaviors of the 
GBM and LAT emission suggests a two component origin 
for the fifth time-bin. Furthermore, the observed spectra of 
GRB 080916C and GRB 090510 show the temporal evolu- 



1/2 

tion rou ghly obeying e P h oc L . from the second ti me-bins 
(see also lGhirlanda et al.|[201flfYonetoku et al.ll2004 ), which 
require a relation 77 oc L 7 ' 16 . For the third and fourth time- 
bins of GRB 090510, when the distinct high-energy compo- 
nent is much brighter than the Band component, we require 
an additional shell behind the two colliding shells and be- 
sides need r a /Y a ~ 10 5 cm close to or even smaller than a 
Schwartzschild radius of mass ~ M@, ~ 3 x 10 5 cm, and 

very small e^es ^3x 10 -6 . 

Other LAT GRBs, such as GRB 080825C (|Abdo et all 

l2009al ) and GRB 081024B jAbdo et al.ll2010h , whose red- 
shifts are not determined, display time-binned spectra fitted 
by Band functions with a temporal behavior of the high- 
energy spectral indices /3 similar to those of GRB 080916C, 
so that they also may be explained by our model. In our 
model the UP emission should exist in all GRBs if the jet is 
in the efficient scattering regime, but for GRBs with small 
L this may be below the detection threshold of LAT, since 
the UP luminosity is proportional to L. This is a general 
statement applicable to wide types of models in which L de- 
termines the flux normalizations of both the Band and the 
high-energy emission components, and consistent with the 
fact that the high-energy emission is detected by L AT only 
in the brightest class of Fermi GRBs (|Granotll2010l ). 

We have shown that if the GRB jet is in the efficient 
scattering regime right from the start, the simple kinematic 
effect causes the onset of the first UP emission to be de- 
layed with respect to that of the first photospheric emis- 
sion by a timescale comparable to the pulse widths or pulse 
separations of the photospheric emission, which should be 
smaller than the variability timescale apparent in the MeV 
light curve (see Section [3]). This may not explain the time 
delays much larger than the apparent variability timescale 
observed in most of the LAT GRBs. In our model, such 
observed large time delays may instead be attributed to a 
temporal evolution of the jet parameters, L, r a /F a , 77, 1Z, 
and t v , a possible origin for which is discussed in the next 
paragraph. From the model fits of the time-binned spectra 
of the three LAT GRBs (in Section [5]), we found that the 
parameter regime should shift from 77 > 77, into 77 < 77, to 
reproduce the observed spectra (this shift leads to the in- 
crease of Z/ up /I/ p h = x; Equation 1181 and 130ft - In this case 
the UP emission component starts with a small delay with 
respect to the photospheric emission onset by the timescale 
comparable to the photospheric variability time, being still 
dim compared with the photospheric emission while 77 > 77, , 
but starts to be bright and detected by LAT when the pa- 
rameters shift into the regime of 77 < 77, . 

This parameter shift appears to be related to the large 
decrease of 7- a /r a , from ~ 10 8 cm to ~ 10 7 cm, in the 
timescale of t d oia y ~ 1 s in long GRBs GRB 080916C 
and 090902B (Table [U). The progenitors of long GRBs are 
thought to be collapsing massive stars, and such a decrease 
of r a /Y a may be explained by the interaction of the jet with 
the stellar env elope before the jet b reakout. The numerical 
simu lations bv lMorsonv et all (|2007l ) (see also lLazzati et all 
2009) show that the jet has three parts around the time 
when it breaks out the star; the jet head, the collimation- 
shocked part, and the free expansion part. The jet head is 
defined as the very thin part at the front end of the jet be- 
tween the forward and reverse shocks, which typically has a 
mildly relativistic speed and the material inside this part es- 
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capes sideways into the cocoon; a large fraction of the jet be- 
hind the jet head, suffering weaker dissipation by collimation 
shocks, has a relativistic speed; and behind the collimation- 
shocked part, the jet material suffers much weaker dissipa- 
tion, expanding adiabatically[f] For the collimation-shocked 
part, r a may be the stellar radius and r a is given by its 
Lorentz factor at the breakout time, while for the free expan- 
sion part, r a may be the size of the central compact object 
and F a ~ 1. In such a jet structure, the emission properties 
may be changed at the transition between the collimation- 
shocked part and the free expansion part. The observed 
time delay would then correspond to the length scale of the 
collimation-shocked portion of the jet at the breakout time, 
Id ~ ctdciay ~ 3 x 10 10 cm. The radius of this part, which 
should be slightly larger than Id, say r a ~ 5 x 10 10 cm, could 
be consistent with the size of the progenitor star and with 
our spectral modeling if F a ~ 500, which is relativistic but 
smaller than i) ~ 3 x 10 3 . The value r a /Y a ~ 10 7 cm for 
the free expansion part is what might be expected from a 
compact object size r a ~ 10 7 cm and F a ~ 1. This r a is a 
few Schwarzschild radii of an object of mass ~ 1OM0. As 
for the short GRB 090510, the parameter shift appears to 
be related to the decrease of r a /T a , from ~ 3 x 10 7 cm to 
~ 3 x 10 6 cm, on a timescale of tdday ~ 0.1 s (Table [2]). 
The model fit indicates an even smaller r a /T a ~ 10 5 cm for 
the third time-bin. This might suggest a smaller size of the 
progenitor and the central object than those of long GRBs. 

In such a scenario, the delay timescale of the LAT 
emission tdelay is related to the global properties of the 
jets and the progenitor systems, rather than the micro- 
physical processes or the simple kinematic (propagation) ef- 
fects, and then it might scale with the total duration of the 
MeV emission td ur - However, in our specific model, while 
t delay ~ f c R*/c scales with the stellar radius R* (where 
f c is the fraction of the collimation-shocked part at the jet 
breakout), tdur may depend on the size of the central region 
of the star R c from whic h the accretion rate is very high, 
i.e., tdur ~ \/ Rc I ( GM C ), where M c is the mass enclosed 
within R c (cf. |K umar et alj 120081 ). Then tdur may signifi- 
cantly depend on the angular velocity and density profiles 
of the star, which could not simply scale with tdelay Ob- 
servations of LAT GRBs with redshifts determined show 
that GRB 080916C, G RB 090902B, and GRB 090926A 
l|Ackermann et all l201ll ) have the durations tdur ~ 15 s, 
~ 8 s, and ~ 4 s, respectively, while the delay timescales 
are all tdelay ~ 1 s. In our scenario tdday is expected to be 
distributed more widely when we obtain a larger number of 
LAT bursts 

The shift from r\ > 77, into 77 < 77, makes the change 
of the photospheric luminosity from ~ L into ~ L(ji/r] t ) a ^ 3 , 
so that it would decrease if L was roughly constant. The 

5 For example, Figure 2 of lLazzati et a! I l l2009h shows that 
the collimation-shocked part has an averaged Lorentz factor of 
~ 10 2 while the free expansion part which follows has a satu- 
rated Lorentz factor of 77 ~ 400. 

6 GRB 090217A d isplays no delay, i.e., t d elav ~ s (while 
*dur(l + z) ~ 33 s lAckermann et alJliblObh . The earliest LAT 
emission may be just the detections of the high-energy tail part of 
the photospheric emission while the UP emission onset could cor- 
respond to the second group of the LAT detections at > 100 MeV, 
~ 7 s after the burst trigger. 



observed MeV emission luminosities appear not to decrease 
(roughly constant for GRB 080916C and GRB 090510 while 
increase by about a factor of two for GRB 090902B), how- 
ever, so that L should increase compensate the shift of the 
parameter regime of 77/77* . Table [2] shows that L increases 
from the first to second time-bins by factors of ~ 3 and ~ 4 
for GRB 080916C and GRB 090510, respectively, while for 
GRB 090902B, the increase factor may be estimated to be 
~ 2/(0.52) 8/3 ~ 10. This may be another parameter tun- 
ing required in our model for LAT GRBs, but it might be a 
selection effect: GRB jets with non-increasing L could have 
dim UP emission which is not detected by LAT. 

For bursts whose LAT emission stays bright from tdelay 
after the trigger, as seen in the three bursts we have studied 
in detail, the parameter regime is assumed to stay 77 < 77, 
in a large fraction of the total duration in our model. This 
means that the radiation efficiencies in those bursts are not 
so extremely large compared with the photospheric emission 
models in which 77 > 77, is assumed over the total duration. 
Indeed we can roughly estimate the radiation efficiencies of 
the three bursts as e 7 ~ £(L p h + L up )tbin/SLtbm ~ 0.4 
for GRB 0809 16C, ~ 0.2 for GRB 090902B, and ~ 0.3 for 
GRB 090510lH Such moderately high efficiencies of prompt 
emission lead to the external shock emission as bright as 
the prompt emission, which is consistent with the observa- 
tions that all the three bursts have bright long-lived high- 
energy emission. For GRB 090902B, its late-time low-energy 
afterglow suggests the total i sotropic kinetic en ergy after 
the prompt phase ~ 10 54 erg (ICenko et alj l20ich . which is 
much smalle r than estimated from the high- energy afterglow 
~ 10 55 erg (|Kumar fc Barniol Durar]|201Ch . This might be 
attributed to the two-components structure of the jet involv- 
ing a narrow and bright spot at the line of sight surrounde d 
by a wider and less energetic region (|Liu fc Wand |2010| ) . 
Many early X-ray afterglows observed by Swtft show very 
steep decays (F u ~ t~ 3 ), which are not compatible with the 
bright high-energy afterglows of LAT GRBs, not decaying 
so rapidly. This is another issue to be solved, and a simulta- 
neous observation by Fermi/LAT and STire/t/XRT from the 
early times will be very helpful. 

The above summary indicates that the photospheric 
emission models may be viable for the temporal and spec- 
tral properties of the MeV and high-energy emission of the 
Fermi/LAT GRBs as well as the other ordinary GRBs, 
although we need the fine tuning of several parameters 
(see a lso the following discussion). Very recently IPe'er et al.l 
(|2010t) showed detailed simulations of the emission processes 
at a dissipation region out of the photosphere of the jet in 
a model similar to ours, including photospheric emission, 
and concluded that such an emission model can explain 
the spectrum of GRB 090 902B, supporting our rough an- 
alytical model. Ilokal (|2010h explores the emission from the 
photosphere and internal shock of the jet with a very high 
77 ~ 10 4 — 10 6 and argues that the high-energy emission of 
the LAT GRBs may be explained as synchrotron emission 
from the internal shock. 



The radiation efficiency of the jet consisting of multiple shells 
coul d be enhanced by the cross IC scatte ring between the shells 
fe.e.. lGruzinov fc Meszarodl200d;lLll2010h . These effects have not 
been considered in this paper. 
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An outstanding problem in the photospheric emission 
models, including ours is that the low-energy spectrum 
of the photospheric emission, usually assumed to be the 
Rayleigh- Jeans part of the blackbody radiation, a pn = 1, 
is much harder than those of the Band components of 
typi cal observed GRBs widely distributed around a ~ 
-1 llPreece et al.ll2000l ; IGhirlanda et alj|2002l ; iKaneko et al.l 
l2006h as well as those of LAT GRBs shown in Figure H 
[6] and [7] In this paper we have mainly considered the two- 
shell system like in Figure [T] but the typical duration for the 
spectral analysis includes many pulses (Note that our model 
parameters and the analysis time-bins are t v (l + z) < 0.04 s 
and T bin > 3 s for GRB 080916C, t v (l + z) ~ 0.3 s and 
T bin ~ 5 s for GRB 090902B, and t v {l + z)~4x 10" 3 s and 
Tbin — 0.1 s for GRB 090510), and thus the superposition 
of the photospheric emission from the multiple shells with 
different e pn has the potential of reproducing the observed 
low-energy spectrum a <C 10 For example, the photospheric 
emissions from multiple shells with different rj and similar 
L, r a /r a , and TZ have a relation L pn oc e pn (not the con- 
ventional blackbody relation L pn oc e ph ; see Equation [4j|, so 
that their superposition appears as a spectrum eF s oc e , 
i.e., a ~ — ljj The parameter rj should be distributed over 
a range of a factor ~ 5 to reproduce the low-energy spec- 
trum extending ~ 1 — 2 decades of energies (jGhirlanda et al.l 
l2007h . However this is just a rough speculation, and we need 
more detailed consistency checks within the model and be- 
tween the model and observations (see also iMizuta et al.l 
l2010t ). The spectral analysis of GRBs with T bin « 0.1 s 
could reveal the intrinsic hard spectrum. Our speculation is 
not inconsistent with the recent analysis that the spectrum 
of as small time-bin as Tbi n ~ 0.5 s ~ t v (l + z) in GRB 
090902B becomes close to a blackbody, while the spectra 
of smaller time-bins but Tbj n > 5 s t v (l + z) in GRB 
0809 16C keep non-thermal (|Zhang et all l20ld ; iRvde et al] 
l2010h . 

The synchrotron (and SSC) emission models also have 
a problem with the low-energy spectral index of the Band 
component. The electrons are required to be cooled so fast 
that the low-energy spectral slope s hould be a = —3/2 i n 
simple types of these models (e.g.. [Ghisellini et alj [2000V 
This is usually softer than the observed slopes, so that the 
superposition effect does not work well. There have been 
many models considering the details of the microphysics in 
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The spectral excesses from the Band component be- 



8 The idea of the superposition of the multiple shells emission 
is similar to a common interpretation of the flat radio spectrum 
of typical blazars (Jp v oc v s with S < 0). The radio synchrotron 
emission is thought to be highly self-absorbed (F v oc ^ 5//2 ), and 
the emission from th e multiple she l ls can reproduce the observed 
spectral slope (e.g., iKonigll Il98ll ; iGhisellini fc Tavecchiol |2009| ; 
iMarschedExjh . 

9 Our spectral modeling have suggested that a correlation rj oc 

L 7 / 16 is required in order to reproduce the observed correlation 
1/2 

£ p h cc L^ h in the time-binned spectra of GRB 080916C and GRB 
090510 (see Section[5]l, which should apply for the maximum value 
of rj in some range of times. 



low ~ 40 keV in GRB 090902B could be one of the im- 
portant discoveries by Fermi/LAT, although such an ex- 
cess was not detected by other detectors extending to suffi- 
ciently low energies s uch as Be ppoSAX and HET E-2 (e.g., 
IGhirlanda et alj l2007h (but see IRvde et"al]|2006h . If these 
are real detections, they may be explained in our model 
as synchrotron or SSC emission from internal shocks, con- 
straining the values of e^es ~ 0.1 (see Table [2J. On the 
other hand, in order for the synchrotron emission not to 
be prominent in the low-energy range in GRB 090510, we 
require e d €s £ 7 x 10~ 3 and < 3 x 10~ 6 for the second 
and third time-bins, respectively. Such a wide range of con- 
strained parameter e_g is similar to the case in the external 
shock, 10 -5 < £b 5- 10~ 2 , constrained from the late after- 
glow observations (jPanaitescu fc Kumar|[2002l \ but it is one 
of the fundamental problems for the shock physics to deduce 
the microphysical parameters es as well as e e and p from the 
first principles. 

The light curves dominated by the synchrotron or SSC 
emission should be correlated with those dominated by the 
UP emission, since they both are produced by the same in- 
ternal shocks. Such a correlation analysis would be useful 
to distinguish our photospheric emissi on model from, e.g., 
other EIC models (|Toma et al.l l2009rj ) . If the high-energy 
emission is produced by up-scattering of soft X-ray emission 
from an external source, such as a cocoon ejected from a 
progenitor star, off the electrons in the jet, the light curve 
of the seed soft X-rays should not be correlated with that 
of the high-energy component. We note that the kinematic 
arguments in Section[3]indicate that a fraction of the photo- 
spheric emission should also be temporally correlated with 
the synchrotron, SSC, and UP emission in our model. 

The synchrotron component could also be detected in 
the optical band. For the model parameters for the second 
time-bin of GRB 090902B shown in Tabled the synchrotron 
self-absorption energy in the observer frame is estimated as 
£syn,a/(l + z) ~ 7 eV and then we have a flux ~ 320 mjy at 
3xl0 14 Hz, which is a very bright optical source even though 
it suffers a strong self-absorption. The bright optical prompt 
emission of some GRBs, being often much brighter than 
the extrapolation of the B and component (e.g. jBriggs et all 
1999; Rac usin et al.ll200ct ). could be attributed to the syn- 
chrotron emission from the internal shocks in our model. 

The spectra of the three LAT GRBs we considered 
may have spectral breaks at e uPt h and at £kn in the high- 
energy range, while a break due to the e ± pair creation in 
the emission site is estimated to be typically much above 
those breaks in our model. Detecting a spectral break in 
the high-energy range by Fermi/LAT and/or by the future 
Cherenkov telescope such as CTA0 would be very helpful 
to constrain the models^ Polarimetric observations of the 
prompt emission also have the potential of distinguishing the 
photosp heric emiss ion models from the synchrotron emission 
models (|Fanll2009l ). In the energy range e < 1 MeV, which 
is the target of some planned missions, the photospheric 



10 http://www.cta-observatory.org. 

11 The distinct high-energy power-law component of GRB 
09092 6A has been co nfirmed to have a spectral break in 
lAckermann et all 112011! ) accepted after the submission of our pa- 
per. 
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emission which consists of the blackbody component and 
the multiple-scattered blackbody component may have very 
low polarization, while some of the synchrotron emission 
models predict de tectable degrees of polarization (see also 
iToma et al.ll2009al . for general discussion of GRB polariza- 
tion). 

The case studies in this paper have shown that the 
photosphere-internal shock model can produce three types 
of the spectral shapes of GRB prompt emission: (i) The 
photospheric emission with the high-energy tail of the pho- 
ton index /3 p h < —2.5 makes the Band spectrum up to 
~ 1 — 10 GeV while a dim UP emission has a small contri- 
bution or just makes the extension of that Band spectrum 
at > 1 - 10 GeV (like the first time-bins of GRB 080916C 
and GRB 090510); (ii) The UP emission is bright but the 
combination of the UP and photospheric components mim- 
ics a Band function with the high-energy index of f3 > /3ph 
(like the second to fifth time-bins of GRB 080916C); (iii) 
The UP emission is bright and makes a high-energy compo- 
nent distinct from a photospheric emission (like the second 
time-bin of GRB 090902B and the second and third time- 
bins of GRB 090510). The increasing number of bursts be- 
ing observed in the LAT field of view (including bursts with 
no detections in the LAT energy range) will determine the 
number ratios of the Band-only spectra and the Band plus 
distinct high-energy component spectra. Extensive paramet- 
ric studies would then be able to clarify whether our model 
with reasonably large parameter space is consistent with the 
high-energy parts of general GRB spectra or only consistent 
with limited numbers of them. A very recent analysis of 52 
bright Fermi/GBM bursts (~ 8 keV - 38 MeV) show that 
the number of the bursts that can be fitted by a Band func- 
tion with 8 > -2.4 is 14 (~ 2 7%) and with B > -2.3 is 7 
(~ 13%) (|Bissaldi et alj|201lh . This implies that the frac- 
tion of bursts with the spectral type (ii) is small, and thus 
we might not require significant fine tuning of parameters 
for general GRBsH 

We have focused on baryonic jets which evolve roughly 
adiabatically and have applied a simple analytical formu- 
lation of the radiation from the photosphere and internal 
shock to deduce the physical parameters of the jets for 
each time-bin. The fine tunings of several parameters ap- 
pear to be required for the three bursts we have studied, 
which might suggest other possibi lities, e.g., non-adiabat ic 
(i.e., significantly dissipative) jets l|Rees fc Meszaros 2005 
or magnetically-dominated jets fe.g.. | Thompsonl 



Spruit et al.ll200ll;lLvutikovll2006l : [Gianniosll2006l ; lFanl 




Zhang fc Yanl 120111 ) Even in those other jet models, how- 



ever, the photospheric emission, its up-scattering effect in 
the internal shocks at large radii^f] and the change of the 
radiation properties due to the interaction of the jet with 
the dense environment could be important for some cases. 



12 The analysis results of the BATSE data (~ 30 keV - 2 MeV) 
show that many spectra are adequately fit ted with a cutoff po wer- 
law function as well as a Band function llKaneko et al.|[2006l) . so 
that they are not suitable for constraining the values and our 
model. 

13 Even in the jets in which the magnetic field energy is dominant 
at the base of the jet, the field energy can be converted into the 
kinetic energy through the jet expansion, leadi ng to the internal 
shocks at large radius (e.g.. iGranot et al.ll201fh . 



Thus, the time-binned spectra and delay timescales of the 
LAT emission onsets could have the potential of revealing 
the radial structure of GRB jets. The difference in the delay 
timescales of long and short GRBs may provide clues for 
understanding the differences between the jets themselves 
as well as their progenitors. 



ACKNOWLEDGMENTS 

We thank the Fermi LAT/GBM group members, 
B. B. Zhang, and T. Sakamoto for useful discussions. 
We also thank the anonymous referee for comments that 
significantly improved the paper. We acknowledge NASA 
NNX09AT72G, NASA NNX08AL40G, and NSF PHY- 
0757155 for partial support. XFW was supported by the 
National Natural Science Foundation of China (grants 
10633040 and 10921063), National Basic Research Program 
of China (973 Program 2009CB824800), and the Special 
Foundation for the Authors of National Excellent Doctorial 
Dissertations of P. R. China by Chinese Academy of 
Sciences. PM is grateful for the hospitality of the Institute 
for Advanced Study, Princeton, during part of this project. 



REFERENCES 

Abdo, A. A., et al. 2009a, ApJ, 707, 580 

Abdo, A. A., et al. 2009b, Science, 323, 1688 

Abdo, A. A., et al. 2009c, ApJ, 706, L138 

Abdo, A. A., et al. 2010, ApJ, 712, 558 

Abramowicz, M. A., Novikov, I. D., & Paczynski, B. 1991, 

ApJ, 369, 175 
Ackermann, M., et al. 2010a, ApJ, 716, 1178 
Ackermann, M., et al. 2010b, ApJ, 717, L127 
Ackermann, M., et al. 2011, ApJ in press (|arXiv:1101.2082[> 
Asano, K., Guiriec, S., & Meszaros, P. 2009a, ApJ, 705, 

L191 

Asano, K., Inoue, S., & Meszaros, P. 2009b, ApJ, 699, 953 

Asano, K., & Terasawa, T. 2009, ApJ, 705, 1714 

Band, D. L., et al. 1993, ApJ, 413, 281 

Beloborodov, A. M. 2005, ApJ, 618, L13 

Beloborodov, A. M. 2010 MNRAS, 407, 1033 

Bissaldi, E., et al. 2011, submitted (|arXiv:1101.3325P 

Bottcher, M., & Dermer, C. D. 1998, ApJ, 499, L131 

Briggs, M. S., et al. 1999, ApJ, 524, 82 

Brunetti, G. 2000, Astropart. Phys., 13, 107 

Cenko, S. B., et al. 2010, submitted [|arXiv:1004.2900|> 

Corsi, A., Guetta, D., & Piro, L. 2010, ApJ 720, 1008 

Daigne, F., Bosnjak, Z., & Dubus, G. 2011, A&A, 526, 110 

Daigne, F., & Mochkovitch, R. 1998, MNRAS, 296, 275 

Dermer, C, & Schlickeiser, R. 1993, ApJ, 416, 458 

De Pasquale, M., et al. 2010, ApJ, 709, L146 

Derishev, E. V., Kocharovsky, V. V., & 

Kocharovsky, VI. V. 2001, A&A, 372, 1071 
Eichler, D., & Levinson, A. 2000, ApJ, 529, 146 
Fan, Y. Z., Piran, T., Narayan, R., & Wei, D. M. 2008, 

MNRAS, 384, 1483 
Fan, Y. Z. 2009, MNRAS, 397, 1539 
Fan, Y. Z. 2010, MNRAS, 403, 483 

Gao, W. H., Mao, J. R., Xu, D., & Fan, Y. Z. 2009, ApJ, 
706, L33 



Photosphere-Internal Shock Model of GRBs 21 



Ghirlanda, G., Celloti, A., & Ghisellini, G. 2002, A&A, 393, 
409 

Ghirlanda, G., et al. 2005, MNRAS, 360, L45 

Ghirlanda, G., Bosnjak, Z., Ghisellini, G., Tavecchio, F., & 

Firmani, C. 2007, MNRAS, 379, 73 
Ghirlanda, G., Nava, L., & Ghisellini, G. 2010, A&A, 511, 

43 

Ghisellini, G., & Celloti, A. 1999, ApJ, 511, L93 
Ghisellini, G., Celloti, A., & Lazzati, D. 2000, MNRAS, 
313, LI 

Ghisellini, G., & Tavecchio, F. 2009, MNRAS, 397, 985 
Ghisellini, G., Ghirlanda G., Nava, L., & Celloti, A. 2010, 

MNRAS, 403, 926 
Giannios, D. 2006, A&A, 457, 763 
Goodman, J. 1986, ApJ, 308, L47 

Granot, J. 2010, in the proceedings of "The Shock- 
ing Universe - Gamma-Ray Bursts and High Energy 
Shock Phenomena", Venice, September 14-18, 2009 
(|arXiv:1003.2452|) 

Granot, J., Cohen- Tanugi, J., & do Couto e Silva, E. 2008, 
ApJ, 677, 92 

Granot, J., & Guetta, D. 2003, ApJ, 2003, 598, Lll 
Granot, J., Komissarov, S. S., & Spitkovsky, A. 2011, MN- 
RAS, 411, 1323 
Greiner, J., et al. 2009, A&A, 498, 89 
Gruzinov, A., & Meszaros, P. 2000, ApJ, 539, L21 
He, H.-N., Wu, X.-F., Tom a, K., Wang, X.-Y., & 

Meszaros, P. 2010, submitted l|arXiv:1009.4716|l 
Ioka, K. 2010, Prog. Theor. Phys., 124, 667 
Ioka, K., Murase, K., Toma, K., Nagataki, S., & Naka- 

mura, T. 2007, ApJ, 670, L77 
Kaneko, Y., et al. 2006, ApJS, 166, 298 
Kobayashi S., Piran T., & Sari R. 1997, ApJ, 490, 92 
Konigl, A. 1981, ApJ, 243, 700 

Kumar, P., & Barniol Duran, R. 2009, MNRAS, 400, L75 
Kumar, P., & Barniol Duran, R. 2010, MNRAS, 409, 226 
Kumar, P., & McMahon, E. 2008, MNRAS, 384, 33 
Kumar, P., Narayan, R., & Johnson, J. L. 2008, MNRAS, 
388, 1729 

Lazzati, D., & Begelman, M. C. 2010, 725, 1137 

Lazzati, D., Morsony, B., & Begelman, M. C. 2009, ApJ, 

700, L47 
Li, Z. 2010, ApJ, 709, 525 
Lithwick, Y., & Sari, R. 2001, ApJ, 555, 540 
Liu, R.-Y., & Wang, X.-Y. 2010, ApJ in press 

(|arXiv:1009.1289|> 
Lyutikov, M. 2006, New Journal of Physics, 8, 119 
Marscher, A. 2009, in Lecture Notes in Physics 794, The Jet 

Paradigm- From Microquasars to Quasars, ed. T. Belloni 

l|arXiv:0909.2576|l 
Medvedev, M. V. 2000, ApJ, 540, 704 
Meszaros, P. 2006, Rep. Prog. Phys., 69, 2259 
Meszaros, P., Laguna, P., & Rees, M. J. 1993, ApJ, 415, 

181 

Meszaros, P.. & Rees, M. J. 1997, ApJ, 476, 232 
Meszaros, P., & Rees, M. J. 2000, ApJ, 530, 292 
Mizuta, A., Nagataki, S., & Aoi, J. 2010, submitted 

(|arXiv:1006.2440|> 
Morsony, B. J., Lazzati, D., & Begelman, M. C. 2007, ApJ, 

665, 569 

Morsony, B. J., Lazzati, D., & Begelman, M. C. 2010, ApJ, 
723, 267 



Murase, K., Toma, K., Yamazaki, R., Nagataki, S., & 

Ioka, K. 2010a, MNRAS, 402, L54 
Murase, K., Toma, K., Yamazaki, R., & Meszaros, P. 

2010b, submitted (|arXiv:1011.0988|l 
Nakar, E., Piran, T., & Sari, R. 2005, ApJ, 635, 516 
Paczynski, B. 1986, ApJ, 308, L43 
Paczynski, B. 1990, ApJ, 363, 218 
Panaitescu A., & Kumar P. 2002, ApJ, 571, 779 
Panaitescu, A., & Meszaros, P. 2000, ApJ, 544, L17 
Pe'er, A., Meszaros, P., & Rees, M. J. 2005, ApJ, 635, 476 
Pe'er, A., Zhang, B. B., Ryde, F., et al. 2010, submitted 

(|arXiv:1007.2228|> 
Pe'er, A., & Waxman, E. 2004, ApJ, 603, LI 
Pe'er, A., & Zhang, B. 2006, ApJ, 653, 454 
Piran, T. 2004, Rev. Mod. Phys., 76, 1143 
Preece, R. D., et al. 2000, ApJS, 126, 19 
Racusin, J. L., et al. 2008, Nature, 455, 183 
Razzaque, S., Dermer, C, & Finke, J. D. 2010, Open As- 
tronomy Journal, 3, 150 
Rees, M. J., & Meszaros, P. 1994, ApJ, 430, L93 
Rees, M. J., & Meszaros, P. 2005, ApJ, 628, 847 
Rosswog, S. 2007, MNRAS, 376, L48 
Ryde, F., et al. 2006, ApJ, 652, 1400 
Ryde, F., et al. 2010, ApJ, 709, L172 
Sari, R., & Esin, A. A. 2001, ApJ, 548, 787 
Sekiguchi, Y., & Shibata, M. 2010, submitted 

[|arXiv:1009.5303|> 
Shemi, A., & Piran, T. 1990, ApJ, 365, L55 
Sikola, M., Begelman, M. C, & Rees, M. J. 1994, ApJ, 421, 
153 

Spruit, H. C, Daigne, F., & Drenkhahn, G., 2001, A&A, 
369, 694 

Thompson, C. 1994, MNRAS, 270, 480 
Thompson, C, Meszaros, P., & Rees, M. J. 2007, ApJ, 666, 
1012 

Toma, K., Sakamoto, T., Zhang, B., et al. 2009a, ApJ, 698, 
1042 

Toma, K., Wu, X. F., & Meszaros, P. 2009b, ApJ, 707, 1404 
Vietri, M. 1997, Phys. Rev. Lett., 78, 4328 
Wang, X. Y., & Meszaros, P. 2006, ApJ, 643, L95 
Wang, X. Y., Li, Z., Dai, Z. G., & Meszaros, P. 2009, ApJ, 
698, L98 

Yonetoku, D., et al. 2004, ApJ, 609, 935 

Zhang, B. 2007, Chin. J. Astron. Astrophys., 7, 1 

Zhang, B., & Pe'er, A. 2009, ApJ, 700, L65 

Zhang, B., & Yan, H. 2011, ApJ, 720, 90 

Zhang, B. B., et al. 2010, ApJ in press (|arXiv: 1009.3 338 ) 

Zou, Y.-C, Fan, Y.-Z., & Piran, T. 2011, ApJ, 726, L2 

This paper has been typeset from a TgX/ P/TgX file prepared 
by the author. 



